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SECTION  1 

TECHNICAL  REPORT  SUMMARY 


1.1  GENERAL 

,m.  roIhe  objective  of  this  research  program  was  to  define  the  preliminary 

Jsd  bror^n8,:y8rea’  ?a8e?  °n  aCOU8tiC  holo^hy  principles,  which 
.  a(l1J  .  ®  detecting  ana  imaging  underground  anomalies  associated  with 

the  ?  ^C4ld  te3t  8lte*  System  requirements  based  on  the  geology  of 

«Jte  were  specified  and  a  proposed  field  test  was  designed. 

niquedwerenmadeUdleS  ^  deVel°P  3  holo8raPhlc  weak  signal  enhancement  tech- 

The  potential  benefits  to  be  gained  from  the  development  of  an  under¬ 
ground  viewing  system  are  many.  For  the  first  time,  a  method  would  be 
available  to  look  into  the  ground  and  "see"  targets,  thereby  defining 
such  anomalies  as  tunnels,  mines,  bunkers,  and  other  voids  or  water-filled 
cavities.  It  would  be  a  valuable  tool  for  surveying  the  earth  ahead  of, 
or  prior  to,  excavations  to  determine  potential  hazards,  for  defining  cer¬ 
tain  fractures,  or  crushed  zones  in  the  earth  such  as  those  which  exist 
after  a  nuclear  underground  explosion,  and  for  locating  underground  bunk- 
ers  and  supply  stores. 


below. 


Specific  tasks  to  be  accomplished  during  this  program  are  given 


1.  Select  a  suitable  field  test  site  at  which  a  future 
field  test  experiment  could  be  performed. 

2.  Perform  a  theoretical  analysis  to  determine  system  require¬ 
ments  based  on  the  geology  of  the  selected  test  site. 

3.  Design  a  preliminary  underground  viewing  system  based 
on  the  system  requirements  determined  from  2  above. 

4.  Perform  an  evaluation  of  presently  available  seismic 
energy  sources  and  detectors  to  determine  if  their 
characteristics  are  consistent  with  system  require¬ 
ments.  Recommend  modifications  of  present  systems, 
or  new  designs,  if  available  equipment  character¬ 
istics  are  not  adequate. 

5.  Investigate  holographic  reconstruction  techniques 
appropriate  for  the  display  of  seismic  holograms, 
with  emphasis  on  the  development  of  image  enhance¬ 
ment  techniques. 

The  results  of  these  tasks  are  summarized  "he  following  sub¬ 
sections.  Detailed  discussions  are  presented  -  ’quent  sections  of 

this  report. 
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1.2  SELECTION  OF  FIELD  TEST  SITE 

The  selection  cf  a  specific  field  test  site  is  a  deviation  of  the 
original  program  plan.  Initially,  Bendix  scientists  and  the  Contracting 
Agent's  Project  Engineer  were  to  jointly  identify  the  types  of  targets 
(tunnels,  miners,  bunkers,  etc.)  of  greatest  interest  to  the  government. 
The  types  of  hard  rock  in  which  the  targets  might  exist  along  with  typi¬ 
cal  target  depths  were  also  to  be  specified.  At  the  first  meeting  be¬ 
tween  Bendix  personnel  and  the  Contracting  Agent's  Project  Engineer, 

Mr.  Gerald  L.  Fitzpatrick,  on  18  June  1971,  it  became  apparent  that  a 
more  practical  approach  would  be  to  identify  the  specific  field  test  site 
at  which  known  and  well-defined  targets  exist  and  where  the  geologic  con¬ 
ditions  of  the  area  were  known  and  documented. 

With  this  modification  of  the  program  plan,  a  search  for  a  potential 
site  was  begun.  Some  of  the  criteria  used  for  site  selection  included 
the  degree  of  surface  weathering,  the  type  of  earth  material  at  the  site 
("hard"  or  crystalline  rock  was  the  principal  geologic  medium  of  interest 
as  specified  in  the  contract),  targets  (tunnels,  caves,  bunkers,  etc.) 
associated  with  the  site,  surface  topography,  and  access  to  the  area. 
Potential  sites  that  were  evaluated  included: 

•  Climax  Stock,  U.S.A.E.C.  Nevada  Test  Range,  Nye  County,  Nevada 

•  Batholith  of  Southern  California 

•  Pikes  Peak  area,  Colorado 

•  Guadalupe  Mountains  area,  New  Mexico  and  Texas 

•  Victorio  Peak  area,  Texas 

•  Big  Bend  area,  Texas  and  Mexico 

•  Winnfield  Salt  Dome,  Winn  Parish,  Louisiana 

A  description  of  these  sites  was  given  in  the  Semiannual  Technical  Report, 
and  is  also  included  in  Section  2  of  this  report  for  completeness. 

The  site  selected  was  the  "Piledriver"  site  in  the  Climax  Stock  on 
the  U.S.A.E.C.  test  range  in  Nevada.  Tentative  selection  was  made  in  a 
meeting  on  September  13,  1971,  with  Bendix  scientists  and  consultants  and 
the  Contracting  Agent's  Project  Engineer.  The  site  was  then  visited  and 
inspected  on  September  14,  1971.  During  this  visit,  A.E.C.  personnel  in¬ 
dicated  that  there  would  be  no  problem  in  arranging  to  use  the  site  for  a 
future  experiment.  They  also  indicated  they  are  very  interested  in  this 
project,  since  it  may  provide  a  means  to  resolve  fracture  zones  caused  by 
nuclear  explosions. 

1.3  ANALYSIS  OF  FIELD  TEST  SITE  TO  DETERMINE  SYSTEM  REQUIREMENTS 

The  geology  of  the  "Piledriver"  site  is  well  documented.  The  rock 
material  consists  of  two  petrological  types,  granodiorlte  and  quartz  mon- 
zonite.  The  two  materials  are  similar  in  both  physical  and  acoustic  pro¬ 
perties.  Velocity  logs  made  at  drill  holes  in  the  area  also  documented 
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this  geology.  Targets  associated  with  the  site  include  a  vertical  shaft 
extending  down  to  two  tunnel  complexes,  one  at  a  depth  of  780  feet,  the 
other  at  1400  feet,  and  chimneys  (fracture  zones)  caused  by  the  nuclear 
explosions.  A  plan  view  of  the  test  site  is  shown  in  Figure  3-1. 

A  model  of  the  test  site  area  was  derived  based  on  known  geological 
information  and  the  results  of  a  refraction  survey  made  in  that  area  in 
1961.  The  model  includes  three  layers  above  a  half-space.  Both  tunnels 
aud  chimneys  formed  by  the  nuclear  explosion  are  located  in  the  half-space. 
The  characteristics  of  these  layers  are  summarized  below. 

•  Near-surface  or  sub-weathering  layer.  Highly  weathered, 
granodiorite,  compressional  wave  velocity  2150  ft/sec, 
approximate  depth  zero  to  20  feet. 

•  Weathered  layer.  Weathered  granodiorite,  compressional  wave 
velocity  3600  ft/sec,  approximate  depth  20  to  100  feet. 

•  Alteration  zone.  Altered  or  fractured  granodiorite, 
compressional  wave  velocity  12,000  to  13,000  ft/sec, 
approximate  depth  100  to  300  feet. 

•  Half-space.  Unaltered  granodiorite,  compressional  wave 
velocity  17,000  to  23,000  ft/sec.  Note:  In  some  cases, 
this  half-space  apparently  starts  fairly  close  to  the 
surface. 

Based  on  this  model  of  the  test  site,  computer  programs  were  written 
to  calculate  ray  paths  in  the  medium,  to  determine  the  maximum  operating 
frequencies  that  can  be  used,  and  to  determine  the  reflection  from  and  mode 
conversions  at  the  various  interfaces.  The  results  from  these  calculations 
were  then  used  in  the  design  of  the  seismic  holographic  experiment. 

1.4  PRELIMINARY  DESIGN  OF  UNDERGROUND  VIEWING  SYSTEM  - 

SEISMIC  HOLOGRAPHY  EXPERIMENT 

The  results  of  the  analysis  of  the  field  test  site  were  used  to 
specify  the  equipment  required  for  a  seismic  holographic  experiment  and 
to  defir.s  the  types  of  seismic  energy  excitation,  excitation  frequencies, 
location  of  the  energy  sources,  and  location  of  and  spacing  of  detectors 
(geophones)  in  the  holographic  array.  Several  constraints  were  adhered 
to  in  this  design,  most  of  which  were  in  the  interest  of  economy.  These 
constraints  included 

•  The  use  of  commercially  available  equipment  (seismic 
energy  sources,  geophones,  amplifiers,  and  recorders), 
if  at  all  possible. 

•  The  use  of  a  reasonably  sized  holographic  aperture  in 
the  vicinity  of  the  surface  over  the  tunnels  where  the 
topography  of  the  area  is  not  extreme. 
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®  The  array  geophones  should  be  located  at  the  surface 
Instead  of  being  buried  in  drill  holes  below  the 
weathered  layer  (approximately  100  feet  deep). 

•  Vibrator  types  of  energy  sources  should  be  operated 
at  the  surface  rather  than  below  the  weathered  layer. 

Calculations  based  on  the  model  described  in  Section  1.3  were  made 
to  predict  the  signal  amplitude  at  the  holographic  aperture  due  to  re¬ 
flections  from  the  target  tunnels  when  ensonified  by  a  monofrequency 
hydraulic  vibrator  seismic  source  at  the  surface.  These  calculations 
were  done  for  excitation  frequencies  of  100,  200,  and  500  Hz.  Calcula¬ 
tions  were  also  made  to  predict  tne  amplitude  of  interfering  reflecting 
and  refracting  waves,  and  an  estimate  of  the  Rayleigh  wave  signal  was 
made.  The  amplitude  of  these  interfering  waves,  compared  to  the  desired 
signal,  determines  the  dynamic  range  aquirements  of  the  recording  equip¬ 
ment.  For  an  excitation  frequency  of  200  Hz  and  below,  the  dynamic  range 
requirements  were  found  to  be  consistent  with  the  available  digital  re¬ 
corders.  Frequencies  above  200  Hz  may  also  be  used,  especially  to  detect 
and  image  the  shallow  (780  feet)  tunnel,  but  a  decision  on  this  should  be 
postponed  until  preliminary  experiments  at  the  test  site  are  made. 

The  spatial  frequency  due  to  signal  reflections  from  the  two  tunnels 
were  calculated.  With  200  Hz  excitation,  reflections  from  the  deep  tunnel 
result  in  a  spatial  frequency  giving  slightly  more  than  four  fringes  across 
a  1500-foot  square  aperture.  (See  Figure  4-1  for  oosition  of  aperture  with 
respect  to  target  locations.)  This  will  result  in  a  rather  poor  resolution 
of  the  reconstructed  image,  but  the  fact  that  the  tunnel  is  long  should 
make  the  image  recognizable.  At  200  Hz,  a  spatial  frequency  giving  approxi¬ 
mately  seven  fringes  will  result  from  reflections  from  the  shallow  tunnel 
(760  feet) .  Here  the  resolution  should  improve  and  a  reasonable  reconstruc¬ 
ted  image  is  expected.  The  use  of  300  Hz  will  give  approximately  10  fringes 
due  to  reflections  from  the  shallow  tunnel,  and  resolution  will  be  further 
Increased.  It  appears  that  300  Hz  is  the  maximum  frequency  that  can  be  used 
with  vibratory-type  energy  sources,  even  for  the  shallow  tunnel. 

The  ability  to  predict  the  amount  of  energy  reflected  from  and  scat¬ 
tered  by  the  chimney  (crushed  zone)  above  the  nuclear  explosion  is  somewhat 
more  difficult.  Although  the  rock  in  the  chimney  area  was  crushed  during 
the  explosion,  and  is  expected  to  be  less  consolidated,  the  velocity  con¬ 
trast,  and  hence  the  effective  reflection  coefficients,  between  it  and  the 
host  rock  is  not  known.  Therefore,  rather  than  an  analytical  prediction  of 
signal  strengths  at  the  surface  due  to  reflection  and  scattering  from  the 
chimney,  an  experimental  determination  during  the  proposed  holographic 
field  experiment  is  recommended.  Since  this  target,  being  primarily  a 
crushed  zone  of  rock,  is  likely  to  be  more  of  a  scatterer  than  a  reflector, 
and  since  it  extends  to  within  500  feet  of  the  surface,  ensonif ication  with 
higher  frequencies  than  that  used  for  detecting  the  tunnels  might  be  feas¬ 
ible.  Ensonif ication  with  an  explosive  source,  and  using  the  highest  use¬ 
ful  quasi-monofrequency  component  of  this  radiation,  is  suggested. 
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Although  it  appeals  that  reasonable  resolution  of  the  shallow  target 
can  be  obtained  when  the  target  is  irradiated  with  a  vibratory  source,  the 
use  of  explosives  as  a  source  of  seismic  energy  was  also  considered.  A 
relatively  .small  charge  of  explosives  (a  few  pounds)  will  induce  many  times 
the  seismic  energy  into  the  earth  than  can  be  obtained  with  vibrator  sources. 
Although  this  type  of  source  is  an  impulsive  source,  ;hat  is,  it  radiates 
over  a  broad  frequency  spectrum,  the  energy  contained  in  a  very  narrow 
(quasi-monofrequency)  band  can  be  greater  than  that  obtained  from  a  mono¬ 
frequency  vibrator.  Proper  design  of  the  cavity  into  which  the  explosive 
is  detonated  can  also  enhance  a  particular  frequency  band.  In  this  man¬ 
ner,  the  intensity  of  energy  in  the  high  frequency  end  of  the  spectrum  can 
be  optimized.  For  this  reason,  the  use  of  explosives,  as  well  as  the  use 
of  a  vibrator  energy  source,  is  recomnended  in  acquisition  of  the  seismic 
hologram.  Rather  than  being  detonated  at  the  surface,  however,  the  ex¬ 
plosives  would  be  set  off  in  a  drill  hole  below  the  weathered  layer  (ICO 
feet  deep) .  Detonating  at  this  depth  has  the  advantage  that  reflected  sig¬ 
nals  from  the  weathered  -  altered  layer  interface  are  virtually  eliminated. 
The  use  of  explosives  as  an  energy  source  does  require  filtering  of  the 
received  signal  to  obtain  a  quasi-monofrequency  signal,  but  this  can  be 
accomplished  during  the  data  processing  in  a  digital  computer. 

To  summarize,  it  is  recommended  that  a  vibratory  seismic  source  be 
used  at  the  surface  to  generate  monofrequency  continuous  and  pulsed  signals 
of  200  H.?  or  greater  for  enaonification  of  various  targets  at  the  site, 
and  that  impulsive  sources  (explosives)  also  be  used  for  target  ensonifica- 
tion.  Further,  it  is  recommended  that  frenuency  modulated  (chirp)  excita¬ 
tion  be  used  to  compare  the  merits  of  these  types  of  excitation  as  related 
to  seismic  holography. 

1.5  EVALUATION  OF  oElSMIC  ENERGY  SOURCES  AND  DETECTORS 

1.5.1  Seismic  Energy  Sources 

Hydraulically  driven  vibratory  sources  have  been  used  in  the 
seismic  exploration  field  for  several  years.  Since  they  are  primarily  used 
for  the  petroleum  industry  and  the  need  is  to  define  the  earth  8  layering 
structure  at  great  depths,  these  sources  are  generally  designed  to  operate 
at  frequencies  frea  about  5  to  less  than  100  Hz. 

Two  approaches  can  be  taken  to  obtaining  higher  frequency 
sources  capable  of  generating  a  monofrequency  signal:  a  new  type  of  source 
might  be  designed,  or  presently  available  vibrators  might  be  modified  for 
higher  frequency  operation.  Since  the  first  approach  would  be  expected  to 
be  more  costly,  an  effort  was  made  to  see  if  the  second  approach  was 
feasible. 

It  was  determined  that  presently  available  vibratots  could  be 
modified  u  operate  as  high  as  500  Hz.  The  modifications  required  for  a 
Wabco  Vibrator  Model  600  B-Df  presently  used  by  United  Geophysical  Corpora¬ 
tion!  a  subsidiary  of  the  Bcndix  Corporation!  include 
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®  Installation  of  a  manifold  and  high  frequency 
servovalve  in  lieu  of  existing  valve 

•  Increase  in  the  ram  area  to  9  square  inches 

•  Installation  of  an  electronic  amplifier  to  drive 
the  new  servovalve 

•  Implementing  a  second  feedback  loop  based  on 
slave  valve  position 

•  Modification  of  ground  contact  pad  to  enhance 
coupling  of  higher  frequencies  into  the  ground 

•  Reinforcement  of  vibrator  structure 

^Ich  these  modifications,  the  displacement  of  the  reaction¬ 
ary  mass  (maximum  displacement  coupled  to  ground)  will  be  5  x  10“5  inches 
at  500  Hz,  2.5  x  10~3  inches  at  200  Hz,  and  HT2  inches  at  100  Hz. 

1*5.2  Seismic  Detectors  -  Geophones 

Commercially  available  geophones  used  in  the  seismic  explora¬ 
tion  field  were  found  to  be  adequate  for  the  holographic  system  require¬ 
ments.  The  minimum  useful  output  of  geophones  is  related  to  the  noise 
level  of  seismic  amplifiers  (typical  amplifier  noise  is  0.1  pV) ;  hence 
geophones  capable  of  detecting  ground  displacements  resulting  in  an 
0.1  pV  output  are  sufficient. 

One  such  geophone  is  the  Geo-Space  Type  GSC-110  detector. 

Its  output  of  0.6  V/in. /sec.  is  constant  over  the  frequency  range  from 
100  to  500  Hz.  This  geophone  will  give  an  output  of  0.1  pV  at  200  Hz  with 
a  ground  displacement  1.35  x  10“10  inch.  Ground  displacements  greater  than 
this  are  expected  due  to  reflections  from  the  deep  tunnel  at  the  selected 
test  site. 

1.6  HOLOGRAPHIC  WEAK  SIGNAL  ENHANCEMENT  STUDIES 

The  development  of  a  holographic  weak  signal  enhancement  technique 
was  initiated  at  Bendix  Research  Laboratories  prior  to  the  award  of  this 
contract.  Further  studies  of  this  technique  as  it  pertains  to  seismic 
holography  were  made  during  this  program.  This  technique,  which  lends  it¬ 
self  to  computer  simulation,  has  shown  promise  in  defining  images  of  tar¬ 
gets  that  have  signal  strengths  too  weak  to  be  seen  by  conventional  holo¬ 
graphic  imaging  (reconstruction)  techniques,  as  is  often  the  case  when  a 
scattering  target  is  located  above  or  below  a  strong  signal-reflecting 
layer. 


Two  computer  models  describing  a  two-layer  representation  of  the 
earth  were  investigated.  The  first  model  contained  a  weak  scattering 
target  below  the  reflecting  interface  of  the  two  layers;  the  second  model 
had  the  weak  target  located  above  the  reflecting  layer.  In  the  computer 


simulations,  all  mode  conversions  (for  example,  P-waves  to  S-waves)  are 
allowed  at  both  the  reflecting  and  scattering  surfaces. 

Results  using  the  first  model  (weak  target  below  the  reflecting 
surface)  show  that  good  image  resolution  can  be  obtained  using  the  en¬ 
hancement  technique  even  when  the  ratio  of  strong  to  weak  signals  is 
1000  to  1.  Using  conventional  reconstruction  techniques,  the  weak  tar¬ 
get  cannot  be  resolved.  Examples  of  these  results  are  included  in 
Section  6. 

Results  with  the  second  model  (weak  target  in  the  zone  above  the 
reflecting  layer)  indicate  that  unattenuated  shear  waves  scattered  from 
the  target  tend  to  interfere  with  the  compressional  waves,  thereby  reduc¬ 
ing  the  resolution  of  the  reconstructed  image.  If  the  shear  waves  are 
allowed  to  attenuate  more  rapidly  than  the  compressionax  waves  (by  a 
factor  of  two  or  greater),  image  resolution  is  improved.  Again,  examples 
of  these  results  are  given  in  Section  6. 

When  the  results  from  the  two  computer  simulation  models  are  com¬ 
pared,  the  image  enhancement  technique  appears  :o  give  better  results 
when  the  weak  target  is  located  below  a  reflecting  plane.  For  th'  nro— 
posed  holographic  field  test,  the  targets  are  located  below  the  rt meet¬ 
ing  layer.  This  will  probably  be  the  case  in  most  field  problems. 


SECTION  2 

SELECTION  OF  FIELD  TEST  SITE 


Th«  original  contract  program  plan  specified  that  Bendix  scientists 
an  the  Contracting  Agent  jointly  identify  the  types  of  targets  (for 

e^m?le'  tunnels»  101068 .  bunkers)  of  greatest  interest  to  the  government 
which  could  be  detected  and  Imaged  using  seismic  holography.  The  type 
of  earth  in  which  these  targets  might  be  located  was  specified  as  "hard" 
or  crystalline  rock.  At  the  first  meeting  with  the  Contracting  Agent's 
Project  Engineer  on  18  June  1971,  it  became  apparent  that  a  far  better 
approach  would  be  to  identify  a  specific  field  test  site  at  which  known 
and  well  defined  targets  existed  and  where  the  geologic  conditions  of 
the  area  were  known  and  well  documented.  With  this  modi  flection  of  the 
program  plan,  a  search  for  a  potential  site  was  begin. 

2.1  CRITERIA  USED  FOR  SITE  SELECTION 

Before  the  search  for  a  specific  field  test  site  was  initiated, 
certain  criteria  were  established  to  be  used  as  guidelines  for  initial 
site  evaluation.  These  guidelines  are  listed  below. 

®  The  area  should  have  little  or  no  soil  overburden  or 

weathering.  This  restraint  will  permit  a  higher  operating 
frequency  to  be  used  for  a  field  test  without  the  need  of 
planting  seismic  detectors  and  sources  below  the  surface. 

•  The  earth  material  at  the  site  should  be  of  a  suitable 
rock  type  (i.e.,  igneous  rock  or  crystalline  limestone) 

to  conform  with  the  contrast  requirement  that  the  principal 
geologic  medium  of  interest  is  hard  or  crystalline  rock. 

•  The  geology  in  u.e  area  of  the  site  should  be  known  and 
reasonably  well  documented. 

•  A  suitable  subsurface  target  must  be  present  at  the  site. 

•  The  area  should  preferably  have  a  reasonably  flat  tonography 
so  that  excessive  phase  compensations  to  correct  for  surface 
elevation  differences  are  not  required  in  an  initial  field 
experiment. 

•  Access  to  the  site  should  be  reasonable  and  there  should  be 
no  restrictions  to  the  use  of  explosives  as  a  seismic  energy 
source . 

2.2  POTENTIAL  SITES  EVALUATED 

Two  Bendix  consultants,  Harry  C.  Kent,  Head  of  the  Department  of 
Geology  at  the  Colorado  School  of  Mines  in  Golden,  Colorado,  and 
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Milton  B.  Dobrin,  Geophysics  Professor  of  the  University  of  Houston 

ln  the 

personrfaalllfr"ltralyh1'’f 

iliar  with  each  site  were  contacted  for  additional  information. 

"Piledriver"6^  W!!k  ^r?f^y  selected  from  the  candidate  sites  was  the 
riledtiver  site  in  the  Climax  Stock  on  the  U.S.A.E.C.  test  ranee  in 

?rE 

personnel  rd'e'onsmian.s'thelon^1^ .  U  ftSSSJ 

was  in^l^dedfiJethI1?tlr  °f  p°tentlal  sltes  that  were  considered 

but  are  also  <  Technlcal  Keport.  Bendix  Report  No.  6050, 

out  are  also  included  here  for  completeness. 

1,2,1  Climax  Stock,  U.S.A.E.C.  Nevada  Test  Rang*. 

Nye  County.  Nevada - 

thp  nnrph  „Th?  Climax  Stock  is  an  intrusive  igneous  body  located  in 
Nevada  Th?»  r«  <P°^  ”  °f,.the  °ak  Sprlng  ')uadriu'81e.  Nye  County, 

is  lmmedlately^north^nf  “c^t™  P°rtl“  °f  ^  TeSt  Slte  °nd 

ihi™  „  _,Th?  araa  0f  exP°sure  of  the  stock  is  approximately  4  square 

ktiLQ°ZTl  ^  uhB  TuCk  tyPeS  are  8ran°diorite  and  quartz  monzo^ite^  The 

?°™Ck  typa?^aVe  been  differentiated  in  field  mapping,  but  there  is 

hv8s° i°  ar89  blHty  dthat  tHe  tW°  r°cks  Wl11  have  essentially  similar 
physical  and  acoustic  properties.  7 

,  uid  a  The  gene,ral  8eologv  of  the  Nevada  Test  Site  is  treated  in 

fP  ^  at;  al,;v  Tha  englueerlng  and  fra^uJa  p^^tlas 

01  the  rock  of  the  Climax  Stock  has  been  studied  by  Ege.2  F 

vertical  shaft^ffe^.lso  ^^Oo^tS:i^^bLdHTe70^^oodt-dOiOle ’t  8 

theUClirted  hemJ8pherlcal  chamber  were  to  be  excavated  in  the  rcclTo" 
the  Climax  stock.  The  vertical  entry  shaft  was  sunk  in  the  proximity  of 

horl^Tj^r®6”  Che  8ranodiorlte  and  the  quartz  monzonlte,  and* the 
horizontal  drift  extended  In  the  quartz  monzonlte  toward  the  m “act 

ho  l“  n“J  drift  on°thelr 

way  :x^.“(r"daeP:al“da;e"8"veritn"s«:iSSt3r  °*  ~* 

H..  Ther,a  ls  aPProxlmately  1000  feet  of  topographic  relief  at 

the  surface  over  the  Climax  stock.  Slopes  ere  steenet  In  tk.  nltlZ 
portion  and  more  gentle  to  the  south  a„dP  so«heL“T«ss  to ‘SfSa 
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should  be  good  because  of  the  numerous  roads  and  trails  put  in  by  the 
Atomic  Energy  Commission. 

2.2.2  Bathollth  of  Southern  California 

Portions  of  Orange,  Riverside,  and  San  Diego  Counties, 
California,  are  underlaid  by  a  complex  of  Igneous  rocks,  mostly  of 
Cretaceous  age,  which  are  known  collectively  as  the  Southern  California 
bathollth.  The  northwestern  portion  of  the  bathollth  has  been  described 
by  Larsen,4  and  other  areas  have  been  reported  upon  by  a  number  of  different 
authors.  Two  areas  within  the  bathollth  seem  to  be  of  prime  Interest: 
the  Cajalco  area  and  the  Cuyamaca  Peak  area. 

Cajalco  Area 

This  potential  site  Is  located  in  Sections  1  and  2,  T. 

4  S . ,  R.  6  W.,  Riverside  County,  California.  In  addition  to  the  publica¬ 
tion  by  Larsen,  a  U.S.  Geological  Survey  open-file  report  by  Page  and 
Thayer6  should  be  available  for  examination  at  the  Menlo  Park,  California 
office  of  the  Geological  Survey.  A  considerable  amount  of  geologic 
mapping  related  to  water  projects  has  also  been  done  for  the  Metropolitan 
Water  District  of  Southern  California,  and  it  should  be  possible  to 
obtain  access  to  these  maps . 

The  principal  rock  of  the  Cajalco  area  is  the  Woodson 

Mountain  Granodiorite.  The  rock  is  typically  white  to  pale  brownish 

gray  with  scattered  black  grains  and  is  rather  coarse  grained.  It 
averages  33  percent  quartz,  60  percent  feldspar,  5  percent  biotite,  and 
small  amounts  of  other  minerals. 

Potential  targets  in  the  area  included  the  old  workings 
of  the  Cajalcc  tin  mine  (this  is  also  referred  to  as  the  Teaescal  mine 
or  Temescal  district  in  some  reports)  and  the  Cajalco  water  tunnel. 

Larsen  (1948)  quotes  older  reports  (esp.  Fairbanks,  1938)  to  the  effect 
that  the  mine  in  the  1890's  had  two  180-foot  vortical  shafts  and  300  feet 

of  horizontal  workings-  The  present  state  of  these  workings.  From  the 

topographic  map  the  Cajalco  tunnel  appears  to  reach  depths  of  100  feet 
or  so  below  the  surface. 

Topographic  relief  in  the  area  is  200  feet  or  less,  and 
access  appears  to  be  good.  Proximity  to  the  Lake  Mathews  reservoir  and 
to  facilities  of  the  water  district  might  restrict  the  use  of  dynamite 
in  testing.  Orchards  now  i ccupy  an  area  adjacent  to  the  tin  mine  pro¬ 
perty,  and  this  might  imply  unfavorable  soil  conditions. 

Cuyamaca  Peak  Area 

The  Cuyamaca  Peak  area  has  been  described  by  Everhart.6 
This  potential  site  is  a  rather  large  area  which  includes  several  more 
restricted  locations  of  interest  within  the  15-minute  Cuyamaca  Peak 
quadrangle.  The  more  important  localities  are  the  Boulder  Creek  district, 
the  Stonewall  Mine,  the  Oriflamme  Mine,  and  .ie  Descanso  Mine. 


Inrated  a  The  Boulder  Creek  district  and  the  Stonewall  Mine  are 
ocated  in  areas  where  the  country  rock  is  a  "mixed"  rock  of  Stonew-ll 

“c"it:osas?b^11rTS?hi8t-  1,16  Jullan  Schi8t  is  a"  oLer  “r- 
eranitic^an  y  ?  T^lasslc  a8e.  «hich  has  been  invaded  by  the 

g  tic  igneous  rocks  of  the  Stonewall  Oranodiorite  which  is  possiblv 

UckUo?BhLaf'<.Thl^  rock  "‘S"‘  problems  because"!?^ 

phlc  rock  bv^the  ^ ’  Ut  1  ,cke  -nva9l°n  and  assimilation  of  the  metamor- 
fai£lv  h  igneous  rock  is  thorough,  the  resultant  might  be  a 

of  homogeneous  rock.  This  would  need  further  investigation.  Many 

of  Everha *"  ^  B°U;!der  Creek  district  were  being  worked  at  the  time 
of  the  ™?',e  s  Paper  or  during  earlier  field  work  in  the  1940 ’s.  Several 

some  200  feet  wo!;kln8S  which  extended  through  a  vertical  distance  of 
.  e  eet  and  3  lateral  distance  of  several  hundred  feet.  The  work- 

to8500ffeet  ?t°neWa}1  M12®  extended  through  600  feet  vertically  and  400 
since  189‘  k  !  S?11  Pr°Perty  haa  been  largely  inactive 

press"  iTcbVlfTr;  bUkrather  l8ree  «*“  °‘  -derate  re.'lel  are 
Thp  b,  |j  r  t8  r^ct*  Stonewall  Mine  is  in  an  area  of  low  relief. 

for  th^  6r  Cf6jk  dlstrlet  would  apparently  be  accessible  for  work  and 

it  lies  iithi^  ynamJte*  The  Stonewall  Mine  might  not  be  feasible  because 

aL  Js  near  d  V  ,  Park’  18  8l0n8  the  8hore  of  Cuyamaca  Reservoir, 
and  is  near  developed  recreational  facilities. 

u  fe.  .  0rlflamrae  Mine  is  located  in  Julian  Schist.  There  raav 

be  sufficient  inhomogeni ties  within  this  rock  to  preclude  ite  use  TheL 
were  originaUy  1000  feet  of  underground  forkings,  but  the  ^ni  is  lone 
abandoned,  and  may  be  largely  caved  ia.  Access  to  the  mine  arL 
would  also  be  difficult  due  to  lack  of  developed  roads. 

rork  4 c  l  j  i  .  The  Descanso  Mine  is  located  in  the  Bonsall  Tonalite.  This 
feldnnnr  lght  gray *  medlum'  to  coarse-grained,  and  consists  of  predominant 
_!*|kP?r»  up  to  25  percent  quart*,  biotite  and  hornblende.  This  is 
probably  the  most  favorable  rock  type  in  the  Cuyamaca  Peak  area  for  the 

andPrhed  e?PCr  Jnient *  ^  mine  originally  had  a  230-foot  inclined  shaft 

nd  three  levels  with  about  80  feet  of  workings  on  each  level.  Access 
to  the  mine  area  should  not  be  difficult,  but  the  workings  may  be  caved. 

2.2.3  Pikes  Peak  Area.  Colorado 

shoui,  P08siblys°::  55  r:br:8^^:ui„-xyr:;„ecoiarada 

“"crrotado  Sptin8a  and  tha  crippia  craak  — 

Peak  bafhnHbhf°th  SlJ68  «e  located  ln  granitic-type  rocks  of  the  Pikes 
lubleorpS  J  ?  complex,  but  the  Cripple  Creek  area  has  been  furthe- 
ubjected  to  volcanic  activity  and  mineralization  related  to  the  fnrmiH 

tend  to^e^h  CrJPple  Creek  caldera*  The  rocks  of  the  NORAD  site  thus 
tend  to  be  rather  homogeneous  while  those  at  Cripple  Creek  are  complex. 
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Pikes  Peek  area  ,Te"^  p'SdTSIr'TT  °l  the 

of^the^cncire°Plkes^Peak  b'ch^H  is' on1«ensi«%tvdy 

site  would  be  the  P;^ciPalfdl8advantage  to  utilization  of  the  NORAD 
Significant fa ?  8teeP  surface  topography.  Topography  would  not  be  a 

reSL  oTr  lln/hVrJPPle  Creek  area*  Extensile  study  of  the 
to  ipfpr  -r  u^!108  *n  C^e  CriPP^e  Creek  district  would  be  necessary 
currently1*1*!  *!!  best  8ite  for  experimentation.  Most  of  the  mines  ale 
desired  depth^an^!  0807  °f  8lS°  deGp  and  Probably  out  of  the 

tfnn  a  i  Sh?uld  the  NORAD  site  be  considered  for  further  investiea- 

dJrine  the^l  °l  8eologic’  engineering  and  mining  data  was  ob^fined 

»o"  d8ba  da,Pl  abLn8tonlr?aVatl0n  °f  the  facilities  It 

desirable  to  obtain  accesp  to  *his  information. 

2,2,4  quadalupe  Mountains  Area.  New  Mexico  and  Texas 
I~hU  rPS”Vr"‘^^  SrTS;  Carlsbad 'caverns  | 

del/1"'  Bcca“se  ot  thc  hl8h  volume  of  visitors  and  tha 
“!f  natu"  ot  thf  '•»  formations,  thc  National  Park  Service  mav  be 
JI2e«  '  PCt"U  the  UM  of  exPloeives .  However,  there  are 

izzZnZtirAzlc°T  ™ and  outside  tha  pa*  -- 

“  -  for^vehicles  ‘“i?™ 

!!'f:'!P  ?”*  15  “'C“SlbU  '"'“-sh  some  of  the  canyons,  and  Le^he 

frot^the  Xai^rr  ‘9  r“88Cd-  *  •*  ^ 


2.2,5  Victoiio  Peak  Area,  Texas 


The  Victorio  Peak  area  is  located  in  westernmost  Texas 
in  Hudspeth  and  Culberson  Counties  and  south  of  the  Cuadalupe  Mountains. 
The  larger  region  around  Victorio  Peak  is  known  as  Sierra  Diablo  and  is 
described  by  King.7 

In  the  vicinity  of  Victorio  Peak  there  is  as  much  as  1500 
to  2000  feet  of  Permian  limestone  of  the  Bone  Spring  and  Victorio  Peak 
formations.  Cave  systems  in  these  limestones  are  not  well  documented, 
but  there  are  vague  indications  on  the  maps  of  the  presence  of  some 
caves.  Some  mining  has  also  been  done  in  the  vicinity  of  small  igneous 
intrusions  which  cut  the  limestones. 

The  principal  drawback  to  the  Victorio  Peak  area  is  the 
topography.  Along  the  east-facing  slope  of  the  Sierra  Diablo  there  is  a 
steep  escarpment.  The  limestones  also  form  cliffs.  The  topographic  and 
geologic  maps  do  show  areas  of  relatively  moderate  relief  along  the  crest 
of  the  Sierra  Diablo  escarpment,  and  if  one  of  these  areas  could  be 
located  in  conjunction  with  a  cave,  the  site  might  be  suitable.  There 
may  also  be  some  difficulty  with  vehicular  access  because  the  road  sys¬ 
tem  is  not  well  developed. 

2.2.6  Big  Bend  Area,  Texas  and  Mexico 

Within  Big  Bend  National  Park  in  west  Texas  and  in  adjacent 
portions  of  Mexico,  thick  limestone  of  Early  Cretaceous  age  are  developed. 
The  geology  of  the  region  is  well  described  by  Maxwell,  et.  al.0  The 
best  areas  of  exposure  ure  in  the  Sierra  del  Carmen  and  its  extensions 
on  the  eastern  side  of  the  Park  (Boquillas  Canyon  area)  and  the  Mesa  <Ie 
Anguila  along  the  western  side  of  the  Park  (Santa  Elena  Canyon  area). 
Similar  rocks  are  also  exposed  in  the  Mariscal  Mountains  in  the  south- 
central  Park. 


Cave  systems  do  exist  within  the  limestones  and  there  have 
been  raining  activities  in  ine  Mariscal  Mountains.  Information  on  the 
cave  systems  is  probably  available  in  speleological  references. 

The  topography  of  the  region  is  rugged,  but  the  tops  of  the 
limestone  mesas  are  relatively  flat.  Vehicular  access  may  be  difficult, 
but  there  are  some  primitive  roads  in  the  area.  Access  to  localities 
on  the  Mexican  side  may  well  be  impossible  except  in  the  vicinity  of 
Boquillas,  Mexico  which  does  have  a  road  connection  to  the  United  States. 

Since  the  areas  are  within  a  National  Park  for  the  most 
part,  it  might  be  difficult  to  obtain  permission  for  the  proposed  experi¬ 
ments.  However,  the  desirable  areas  are  /emote  where  the  visitor  use  is 
low,  so  permission  might  be  granted  with  less  difficulty  than  in  other 
National  Parks.  There  is  a  slight  possibility  that  a  suitable  site  might 
be  located  outside  the  Park  boundary,  but  the  best  locations  appear  to 
be  within  the  Park. 
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2.2.7 


Although  salt  is  not  a  hard  rock  material,  lts 
properties  (high  velocity)  are  similar.  For  this  reason,  salt-bearing 
areas  were  not  excluded  from  the  survey. 

The  Minefield  salt  dome  Is  located  in  Winn  Perish,  loulslene, 
approximately  midway  between  Shreveport  and  Alexandria.  _.The.salt,  which 
extends  within  300  feet  of  the  surface,  has  been  mined “ * ‘l£?  klngs 
at  a  depth  of  about  800  feet.  Maps  are  available  of  the  mine  working 
The  cap  rock  on  the  salt  has  been  quarried,  having  a  layer  of  gyps 
and  anhydrite  (which  have  physical  properties  similar  to  rock  sal  ) 
between  the  surface  and  the  salt. 

Unfortunately,  the  tunnels  and  mine  workings  have  been 
flooded  with  wate’  in  the  past  few  years  and  the  mine  is  no 
ible.  In  addition,  the  area  quarried  over  the  salt  is  pr  y 
than  would  be  desired  for  a  field  test. 
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SECTION  3 


ANALYSIS  OF  VIEWING  SYSTEM  REQUIREMENTS  BASED 
ON  GEOLOGY  OF  FIELD  TEST  SITE 


In  this  section,  the  geology  of  the  selected  field  test  site  is 
described  in  more  detail  and  a  model  of  the  site  is  formulated.  The 
seismic  velocities  (for  compressional,  shear,  and  surface  waves)  as 
well  as  Poisson's  ratio  are  established  for  all  zones  in  the  model. 

The  attenuation  of  seismic  waves  in  the  different  layers  are  calculated. 
Calculations  giving  the  response  of  the  tunnels  to  seismic  wave  irradia¬ 
tion  are  also  made. 

3.1  TOPOGRAPHIC  AND  GEOLOGIC  FEATURES  OF  THE  CLIMAX  STOCK,  AREA  15, 

NEVADA  TEST  SITE,  NEVADA 

The  primary  geologic  feature  of  Area  15,  the  Climax  Stock,  is 
shown  on  the  enclosed  geologic  map,  1-328,  published  by  the  U.S.  Geo- 
] jgical  Survey.  An  index  map  locating  the  Nevada  Test  Site  is  also 
given . 

Figure  3-1  details  existing  tunnels.  All  tunnels,  entrance  shafts 
and  nuclear  detonation  areas  from  Hardhat  and  Piledriver  experiments 
lie  within  the  Climax  Stock.  This  Stock  is  composed  of  two  distinct 
petrological  types  in  the  immediate  area  of  the  target  tunnels.  The 
lithologic  contact  strikes  irregularily  west  to  northwest,  with  the 
granodiorite  member  to  the  north  and  quartz  monzonite  to  the  south. 

The  granodiorite  is  oldest  with  the  younger  porphyritic  quartz  monzonite 
aged  230  +  25  million  years.  Map  1-328,  sheet  2  of  2,  contains  sections 
through  the  Climax  Stock.  The  Stock  is  thought  to  extend  beyond  13000 
feet  in  depth  with  the  granodiorite,  quartz  monzonite  contact  dipping 
steeply  to  the  south. 

Drill  holes  and  coring  in  the  area  indicate  considerable  jointing. 
No  core  recovery  from  depths  above  20  feet  is  reported  from  any  of  the 
core  holes.  During  the  inspection  of  this  site  on  14  September  1971, 
core  samples  from  a  drill  hole  were  viewed  which  indicated  competent 
rock  within  10  feet  of  the  surface.  Apparently,  the  combination  of 
fracturing  and  weathering  of  rocks  near  the  surface  and  drilling  consid¬ 
erations  have  combined  to  inhibit  core  recovery  at  these  shallow  depths. 

Some  faulting  or  fracture  zones  exist  as  can  be  seen  from  the  E-E' 
and  F-F'  sections  which  contain  the  U-15a  exploratory  dr_ll  hole. 

Several  joint  sets  have  been  mapped  from  drill  cores.  Detailed^ 
information  can  be  obtained  in  the  Geological  Survey,  TEM  Report  836. ‘ 
Some  jointing  and  fracture  zones  are  extensive  and  can  be  correlated 
between  drill  holes  (See  Figure  3-2  which  is  a  section  through  the  drill 
holes  mapped  in  Figure  3-1.) 
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Figure  3-1  -  Map  Showing  Location  of  Piledriver  Tunnels 
Area  15,  Nevada  Test  Site,  Nevada 

and  Granite  VrtnSll™1™1  S'“  °f  J°lntS  ln  tha  araa  °f  «“ 

(1)  Strike  N30  to  40  degrees  W,  dip  15  to  35  degrees  NE 

(2)  Strike  generally  NW,  steep  dip 

(3)  Strike  generally  NE,  steep  dip 

**  l0W  dlp  a"8le  a11  8eem  to  be  flrn,ly  healed  with  secondary 
minerals  (quartz,  pyrite,  and  feldspar).  The  joints  with  nearly  vertical 
dip  are  less  well  healed,  but  are  generally  filled  with  clay  and  calcite 
Velocity  surveys  run  from  shot  points  1000  feet  away  from  the  U-15a  drill 
hole  show  that  average  horizontal  velocities  are  slightly  higher  than 

vIrMflV:r,tlCal  VeloClbles-  indicating  that  open  oAoosely filled 
tical  joints  are  probably  not  an  important  factor.  The  number  of 
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INDICATES  MAPPED  FRACTURE 
AND/OR  FAULT  ZONE 


Figure  3-2  -  Fault  Zone  Plunge  Section 


joints  decreases  with  depth,  especially  below  about  400  feet  in  the  U-15a 
area,  and  the  dip  of  the  joints  generally  increases  with  depth.  In  fact, 
the  extent  of  in  situ  jointing  and  fracture  zones  decreases  with  depth 
to  a  sufficient  degree  that  detection  of  distressed  fractured  zones 
surrounding  the  tunnels  has  been  reported. ^ ^ 

Additional  geologic  information  consisting  primarily  of  density 
and  porosity  data  and  the  results  of  geophysical  logs  is  given  in  a 
Geological  Survey  report  by  Roach,  Izett,  and  Roller,  June  1959, 12 
Interpretation  of  Geophysical  Logs  of  the  Granite  U-15a,  Dolomite  Hxll 
No.  1  and  Marble  No.  Me-2  Drill  Holes,  Nevada  Test  Site,  Nye  County, 
Nevada.  The  summary  of  the  report  indicates  negligible  variation  in 
density  or  porosity  between  the  two  rock  types.  Unfortunately,  data  on 
the  seisTiically  important  porosity  are  insufficient  closer  than  144 
feet  from  the  surface;  however,  the  percent  porosity  of  weathered  surface 
samples  of  both  the  quartz  monzonite  and  the  granodiorite  is  higher  than 
that  of  the  unweathered  rock.  The  grain  density  of  both  weathered  and 
unweathered  material  is  essentially  the  same,  and  the  dry  bulk  density 
of  both  is  nearly  the  same,  although  weathered  values  of  dry  bulk  density 
appear  to  average  less  than  those  for  unweathered  material. 

Weathering  effects  can  extend  to  as  much  as  360  feet  in  the  drill 
holes  studied.  In  the  upper  portions  of  the  drill  holes,  there  is  both 
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"Granite ""hole 'to  a^^t.1"111"318  t0  ^  ^  be  effective  in  tha 

in  fho^6  C°mplex  of  tunnels  in  the  Climax  Stock  lies  almost  coraplt'elv 
in  the  quartz  monzonite  portion.  The  Hardhat  drift  lies  at  a  depth  oJ 

shaft^lSOO  shaft) iwill1Ip^vtUnnel  SyStem  3£  14°°  feer‘  ^  entrance 
t-ho  f  '  ^  serve  as  a  vertically  oriented  target  as  will 

the  fracture  zone  surrounding  the  Flledrlver  detonation  aria' 

3.2  COMPILATION  OF  GEOPHYSICAL  INFORMATION  CLIMAX  STOCK.  AREA  15 

Area  15I’are°reonrJj!eiS','1C  Inv“ “nations  which  have  been  carried  out  at 
Logan- Lolllpop^°Project  ??“  iLllTtlTM 

attenuation  measurements,  and  some  uphole  velocity  plots. 

centra i!8U!ie  lndl‘:ate8  the  area  covered  by  these  surveys.  They  lie 

Kmni  *  f  tHe  "°rth  °f  the  Hardhat  aad  Piledriver  tunnels  The 

Piledrive^detonatlon^06  1  fr*""  3'4)  -  •*»  ""  the 

inform^M  Velk  Un±f°rm  rePorts  have  provided  a  substantial  amount  of 
tha  ^  if"  abrUt  thG  se±smlc  Properties  of  the  Climax  Stock;  however 
an  extranilat?  L°llipop  survey  lines  and  drill  hole  data  requires 

grap^  exDerWntrUtiWa  proPosed  site  *<>r  the  seismic  holo- 

ful-  in  narM  ?  <  saarch  for  more  inforniation  has  not  been  success- 

sii«  Jt  is  a  ?ocI;  "r  ",'?  fro"  the  1500  entrance  shaft  are  desirable 

no  velocity  data  wlthin"l50  fee?  of  ^surface.  Al8°’  Chere  *”  ',ll”OS,; 

la»er0"ve^M“e"  ^  15  to  25  feet  ln  Thickness.  SelL  this 

va^I’  Ve-ocl5les  average  approximately  12500  ft/sec,  although  observed 
alues  range  between  approximately  11200  and  13500  ft/sec.  These 

d^^erences  are  undoubtedly  due  to  weatherine  altpi-apinn  j  , 
changes  in  norosim  ,  eacnenng,  alteration,  and  subsequent 

usual! v  Calf?  !  T  y’u  Correlatlon  between  velocity  and  composition  is 
sualiy  valid  on.  y  when  porosity  is  reduced,  i.e.,  at  high  pressure 

~  f3^®"  W*|lch  affect  velocities  in  rock  are  mineral  composition 

»Uh  »ana«  «  StratUrei.Prf8SUr'-  *ral"  slaa>  cementation  “™c?L 

Tn  L?  ?!  bedding  or  foliation,  and  alteration.  However  scatter 

porosity  ^he  intJi  d?pths  anJ  low  Pressure  can  usually  be  ascribed  to 

ss 
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Figure  3-3  -  Map  of  Climax  Stock  Area 

Additional  data  are  obtained  from  uphole  surveys  in  the  area  such 
as  Figure  3-5  which  is  an  uphole  velocity  plot  for  the  L-2  drill  hole, 
the  nearest  to  the  proposed  experiments.  The  entire  hole  lies  within 
the  granodlorite,  and  the  values  can  only  suggest  the  chart.:  ceris  tics 
of  the  quartz  monzonlte.  The  velocity  profile  suggests  a  flnr.r  weather¬ 
ing  structure  than  noted  in  the  refraction  survey.  The  two  zones  will 
be  termed  the  near-surface  layer  (NSL)  with  a  thickness  of  approximately 
20  feet  and  the  weathered  layer  (WL)  with  a  thickness  of  approximately 
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Figure  3-4  -  Lollipop  NSVP  Line  No.  1 

80  feet.  The  high  speed  refraction  12900  ft/sec  would  correspond  to  a 
refraction  in  the  layer  below  (WL) ,  i.e.,  an  altered  layer  (AL)  in 
which  some  lithologic  change  has  taken  place  but  weathering  is  not 
predominant. 

The  geologic  report,  "'EM  Report  836 gives  detailed  petrologic 
information  on  the  U-15a  Exploration  Hole.  This  hole  penetrates  both 
rock  types.  Generally,  quartz  monzonlte  extends  to  325  feet  and  grano- 
dlorlte  was  occasionally  encountered  in  the  0  to  325  feet  zone.  The 
geologic  report  yields  evidence  of  hydrothermal  alteration  of  the  quartz 
monzonlte  and  granodlorlte.  The  alteration  products  Include  clay 
minerals,  chlorite,  secondary  feldspar,  serlclte,  quartz,  epldote  and 
sulfide  minerals,  mainly  iron  sulfide.  The  clay  minerals  and  chlorite 
occur  In  zones  throughout  the  rock  and  are  concet  trated  with  other 
alteration  products  along  prominent  northeast  dipping  Joints.  The  clay 
minerals,  predominantly  ruUiitmorlllonlte ,  constitute  only  a  fraction 
of  the  total  rock  and  appear  to  be  more  common  In  zones  1  to  15  feet 
thick  associated  wiU«  9teeply  dipping  fractures.  The  minerals  formed 
by  hydrothermal  alteration  were  abundantly  distributed  from  the  surface 
to  a  depth  of  about  360  feet.  In  this  zone,  both  the  granodlorlte  and 
quartz  monzonlte  contain  portions  In  which  the  alteration  Is  extensive. 
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Figure  3-5  -  Uphole  Tine  Depth  Curves 


A  summary  of  the  above  information  suggests  the  following  model  of 
the  seismic  characteristics  for  whe  Climax  Stock. 

Layer  1.  Near-Surface  Layer 

Thickness  20  feet  +  10  feet.  Composition  highly  weathered  quartz 
monzonlte  and  alluvium.  The  alluvium  referred  to  here  is  the  result 
of  several  low  southeast-trending  ridges  which  have  40  to  60  feet  of 
relief.  The  intervening  valleys  contain  Intermittent  streams  and  are 
partly  filled  with  rock  debris  derived  by  weathering  of  the  Climax 
intrusive  mass.  Velocity  of  compresslonal  waves  is  2150  ft/sec  and 
density  is  5.15  slug/ft3.  The  porosity  is  undetermined  but  certainly 
greater  than  3  percent. 
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Layer  2.  Weathered  Layer 

Thickness  80  feet  +  20  feet.  Composition,  weathered  and  hydrother- 
mally  altered  quartz  oonzonlte.  In  this  layer,  the  weathering  effect 
upon  elastic  wave  propagation  is  predominant.  The  quartz  monzonlte  is 
medium  gray  to  light  gray  and  contains  phenocrysts  in  a  fine-to-medium 
grained  ground  mass .  Occasional  penetrations  of  granodlorlte  may  occur 
in  this  zone,  particularly  over  the  northernmost  end  of  the  tunnels. 

The  granodlorlte  Is  mostly  light  gray  to  greenish  gray,  equigranular ,  and 
medium  grained.  However,  this  difference  is  not  expected  to  change 
velocity  significantly  in  the  two  rock  types,  particularly  when  they 
are  highly  weathered.  Compresslonal  wave  velocity  is  360C  ft/rec. 

Layer  3.  Altered  Layer 

Thickness  200  feet  +  50  feet  •  This  layer  is  characterized  by  a 
significant  compresslonal  wave  velocity  change  when  compared  with  the 
weathered  layer.  Average  velocity  is  12900  ft /sec.  A  significant  velo¬ 
city  scatter  is  observed  from  the  refraction  data,  which  probably 
indicates  variation  in  the  thickness  and  extent  of  alteration  nones. 

This  layer  will  act  as  a  strong  reflector  when  surface  sources  are  used. 

Layer  4.  Unaltered  Layer 

This  layer  is  sufficiently  thick  so  that  it  may  be  considered  a 
half-space.  Data  obtained  from  seismic  studies  in  the  Piledrlver  tunnels 
yield  an  average  velocity  of  19400  ft/sec.  Density  is  5.15  slug/ft^. 

This  implies  no  substantial  difference  in  density  for  the  entire  test  site. 
A  destredsed  zone  consisting  of  fractured  and  jointed  rock  has  been 
observed  surrounding  the  Piledrlver  tunnels.  From  uphole  surveys  and 
refraction  shooting,  a  range  of  values  of  velocity  and  depth  of  the  zone 
has  been  reported.  Velocity  ranges  from  4300  to  15200  ft/sec  and  the 
thickness  of  destressed  zone  from  0.3  to  19.8  feet. 

3.3  EVALUATION  OF  SHEAR  WAVE  VELOCITIES  AND  POSSION'S  RATIO  FOR  ALL 

ZONES  AT  NEVADA  TEST  SITE 

A  tabulation  of  measured  shear  wave  velocities  for  igneous  rocks 
can  be  found  in  Clark,  Handbook  of  Physical  Constants.*4  The  range  of 
values  for  quartz  monzonlte  and  granodlorlte  are  9500  to  10500  ft/sec. 

When  pressure  and  temperature  are  low,  corresponding  compresslonal  wave 
velocities  are  15100  to  17250  ft/sec  which  suggests  rather  competent 
rock,  comparable  to  the  altered  or  unaltered  layers  found  at  the  Nevada 
Test  Site.  The  Poisson's  ratio  calculated  from  the  above  velocities 
lies  in  the  range  0.17  to  0.21.  In  general,  the  value  of  Poisson's  ratio 
will  Increase  for  less  competent  rock,  but  data  from  rock  with  a  greater 
burial  depth  indicate  that  Poisson's  ratio  will  rarely  decrease  below  a 
value  of  0.17.  If  value  of  0.2  is  assumed  for  both  altered  and  unaltered 
layers,  then  the  respective  shear  velocities  will  be  7900  and  11900 
ft/sec. 
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3.4  ATTENUATION 

effnrrD!IIln8  ?*  refrac5lon  experiments13  at  the  Nevada  Test  Site,  an 
usfno  f?  ^“ad!!  t0  ohtaln  information  on  compressional  wave  atcenuation 
S? L “S?  8aln  recordi"8  of  the  24  amplifiers  in  the  geophone  spread 
The  amplifiers  were  also  tapered  so  that  amplifier  sensitivity  was  in- ‘ 

ZlVtLat  8r?ft£!  dl8tancea  ^om  the  source .  From  a  paper  rLoVd  Vhe 

the TtT  mea8"ed  for  the  first  or  seVoVd  excursion  on 
th  in  f'u  j  resulting  amplitude  responses  were  plotted  in  db  versus 

orovlrif  °f  horizontaJ  distance  from  shot  point.  The  slope  of  these  plots 

?raca  aL?<rrUrV!,a"  lnV8r8a  p0“ar  relationship  between  salanogram 

™natS  In  thaaS«f  ?ta"C8  'j0”  tha  8h0t  poln,:  °rl8‘"-  Tha  «ati.cical 
i  ..  ta  is  very  large  and  the  derived  attenuation  values 

can  oniy  provide  an  estimate  of  the  exact  conditions.  Particularly  un¬ 
fortunate  was  the  lack  of  frequency  information  since  attenuation  is 

the^range  'oIo  sTh.  "“fT}  "'"“l0”  aurva>’a  “•  1°»  fra, nancy,  in 
range  zu  to  oU  Hz#  It  is  assumed  in  this  mao  *.u«  j_ 

PimlLV  rdtotal  >a"l"»  . 

i  ..j  *  ,  8ure  3  6  includes  applicable  plots  of  attenuation  from 

rabial-l  fraCtl°n  L1"a  N°-  U  Avara*a  va  uee  «f  slope  ara  glvan  ln 
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Table  3-1  -  Average  Slope  Values 


Spread 

Slope 

Velocity 

Rock  Type 

I,12-NSV13+ 

1.8 

13100  ft/sec 

Granodiorite 

L11-NSV12 

1.5 

13500 

Granodiorite 

L13-NSV12 

2.0 

12800 

Granodiorite 

L13-NSV14 

2.1 

13750 

Granodiorite 

?  * 

10200 

Quartz  monzonite 

L14-NSV13 

1.5 

11200 

Granodiorite  and 

Quartz  monzonite 

L15-NSV14 

2.9 

14500 

Quartz  monzonite 

L12-NSV13  means  shot  point  at  L12,  spread  extends  toward  L13  along 
Lollipop  Line  No.  1.  Shot  point  separation  is  74C  feet. 

The  break  in  "ttenuatlon  on  L13-NSV14  is  probably  due  to  the  change  in 
lithology  from  granodiorite  to  quartz  monzonite.  There  is  also  a 
change  to  a  lower  velocity  along  the  refraction  line. 

From  this  table,  the  average  attenuation  slope  is  1.8  for  grano¬ 
diorite  and  the  average  compresslonal  wave  speed  is  12900  ft/sec.  For 
quartz  monzonite,  the  slope  is  2.2  and  the  average  velocity  12000  ft/sec. 
This  suggests  a  generally  higher  attenuation  and  lower  velocity  in  quartz 
monzonite,  which  is  to  be  expected  when  the  differences  in  petrology  of 
the  two  rock  types  are  taken  into  consideration.  However,  statistically 
the  error  inherent  in  the  quartz  mon:  onlte  data  precludes  a  definite 
statement.  Therefore,  the  average  values  obtained  from  the  granodiorite 
will  be  used  as  representatives  of  the  altered  layer  for  the  entire 
Nevada  Test  Site. 

Attenuation  of  elastic  wave  energy  in  rock  is  known  to  be  distance 
dependent  of  the  form  exp(-uy),  where  a  is  frequency  dependent.  Attewell 
and  Ramana,15  "Wave  Attenuation  and  Internal  Friction  as  Functions  of 
Frequency  in  Rocks"  have  compiled  values  of  a,  and  White16  lists  specific 
experimental  values  in  Seismic  Waves,  Radiation,  Transmission  and 
Attenuation.  Figure  3-6  contains  a  series  of  semi logarithmic  plots  of 
attenuation.  A  straight  line  plot  is  described  by  the  linear  form 
equation 


Y  -  Bx 
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where 


Y  -  log  y 

y  is  horizontal  distance  In  feet 
x  Is  decibels  (dB) 

B  Is  slope 

Changing  decibel  to  nepers  yields 

log  y  •  B  (8.686)  n 


where  n  is  nepers.  Slope  is  -1.8  and 

log  y  •  15.65  (-n) 

Values  for  n  are 


y 

n 

dB 

a  •  dB/y 

500  ft. 

0.1725 

1.5 

0.003 

1000  ft. 

0.192 

1.67 

0.00167 

1500  ft. 

0.203 

1.765 

0.00118 

value  of  a  ■ 

2  x  10"3  dB/ft. 

Following  the  assumption  that  first  arrivals  on  the  record  have  a 
frequency  of  approximately  40  Hz,  then  a  •  a*f  is 

a  -  5  x  10-5  f  ~ 


where  f  is  frequency  in  Hertz. 

This  value  of  attenuation  is  of  the  correct  order  of  magnitude 
when  compared  with  published  values.15*16  and  corresponds  to  the  attenua¬ 
tion  for  the  altered  layer. 

Once  this  value  for  the  altered  layer  has  been  obtained,  rough 
estimates  for  the  attenuation  in  the  other  layers  can  be  made  by  scaling 
the  data  based  upon  the  assumption  that  equal  wavelengths  have  equal 
attenuation.  Such  an  assumption  indirectly  takes  into  account  the 
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important  weathering,  alteration,  and  porosity  changes  which  occur  in 
the  rock.  These  factors  determine  the  wave  speed  and  therefore  also  the 

wavelength. 

First  the  values  of  a  for  100  Hz,  200  Hz  and  500  Hz  in  the  altered 
layer  are  calculated 


For  the 


a100  Hz 
°200  Hz 
°500  Hz 

unaltered  half  space. 


-  5  x  10"3  dB/ft 

-  1  x  lO"7  dB/ft 

-  2.5  x  10"2  dB/ft 
the  scaled  values  are 


a100  Hz 
a200  Hz 
a5u0  Hz 

For  the  near  surface  layer 


3.3  x  10“3  dB/ft 
6.6  x  10"3  dB/ft 
1.65  x  10-2  dB/ft 


a100  Hz 
a200  Hz 
a500  Hz 


3.1  x  10"2  dB/ft 

6.1  x  10'2  dB/ft 
1.525  x  10_1  dB/ft 


For  the  weathered  layer 


a100  Hz 


a200  Hz 


a500  Hz 


1.8  x  10-2  dB/ft 
3.6  x  10'2  dB/ft 
9.0  x  10-2  dB/ft 
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vpinr<fince  ®8timates  are  now  available  for  congressional  velocity,  shear 
velocity,  and  attenuation,  it  is  possible  to  calculate  reflection  and 
transmission  coefficients  for  ell  layers  at  the  Nevada  Usl  sllT. 
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nevadaTtestAsiteRANSMISSI0N  coefficients  for  layering  at  the 


vihratnr  th®  8U^face*  the  source  of  wave  energy  will  be  a  hydraulic 

forces  actip8  on  the  surface  of  the  ground  will  be  oriented 

atteiuates’withVd/rPa8a5i0n  characterlzed  by  a  displacement  which 
attenuates  with  distance  from  the  source.  In  the  near-surface  layer, 

^00  „®veJen8th  of  the  500  Hz  compressional  wave  is  4.03  feet,  and/or  the 
H*  8  *  waves  is  1.72  feet.  These  wavelengths  are  of  the  orde^of 

with  dimpn8j°n8  °f  vibrator  8 r°und  contact  pad.  For  a  contact  pad 
Innut-  Jor  !  !  compared  to  the  wavelength,  the  relationship  between 

input  force  and  output  wave  displacement  is  well  known.16 

For  the  compressional  wave  displaceiuent  u 


G  cos  6 
o 

1  -  2  (V  /V  )2  sin2  & 

■ —  s  p 

-iur/v 
e  ** 

i  iwt 
e 

*">vp2*  ( i 

1  ’  2  <VVn>2 

L  s  P 

:  2  1 

sin  6 

2 

+  4  (Vs/Vp)3  sin2  6  cos  6  [l  -  (Vfl/Vp)2  sin2  s] 


For  the  shear  wave  displacement  u 

$ 

„  s  n  r  2  O  1/2  -iwr/V 

u  ^  ~G0  8in  6  cos  6  |(V8/Vp)  -  sin2  6  ]  e  8  eiwt 

4*  r  .  .  . 

npV8  r^(l  -  2  sin2  6)2  +  4  sin2  6  cos  6  [(Vg/Vp)2  -  sin2  6  ]  ^  } 


where  GQeia,t  i8  the  input  point  force 
r  is  radius 

<5  is  angle  from  vertical 
(i)  is  radians  (source) 

Vp  is  compressional  wave  speed 
Vg  is  shear  wave  speed 
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For  reflection  and  holographic  experiments,  the  wave  energy  entering 
at  almost  vertical  angles  of  incidence  is  of  greatest  importance.  Plot¬ 
ting  the  curves  for  ur  and  u<*,  shows  immediately  that  ur  is  the  predominant 
wave  mode  in  this  range.  In  addition,  reducing  the  pad-to-wavelength 
ratio  results  in  increasing  the  directional  characteristic  so  that  even 
more  compressional  wave  energy  is  concentrated  into  a  narrow  cone  about 
the  vertical.  If  then  the  u^  component  is  neglected,  only  the  ur  com¬ 
ponent  will  arrive  at  the  first  boundary  between  NSL  and  WL.  The  ur 
displacement  can  be  divided  into  horizontal  and  vertical  components. 

Indeed  for  angles  of  incidence  close  to  vertical,  ur  approximates  ux, 
the  plane  wave  displacement. 

Transmission  and  reflection  coefficients  can  be  obtained  for  the 
various  layers,  based  upon  an  analysis  by  Muskat  and  Meres.17  The 
coefficients  are  based  upon  the  propagation  of  plane  waves  in  an  elastic 
medium.  Therefore,  this  theory  only  approximates  the  wave  front  obtained 
from  the  vibrator  source.  The  approximation  is  improved  however,  as  the 
passes  through  successive  layers,  an  effect  discussed  by  Cagniard,18 
Reflection  and  Refraction  of  Progressive  Seismic  Waves. 

Whenever  waves  are  incident  upon  an  interface,  mode  conversion  will 
occur.  This  splitting  of  wave  energy  at  every  interface  yields  a  pro¬ 
liferation  of  waves  which  are  a  major  complication  in  record  interpreta¬ 
tion.  At  steep  entrance  angles,  the  compressional  wave  dominates  where¬ 
as  at  angles  away  from  the  vertical,  downgoing  shear  waves  from  the 
surface  source  become  important  and  at  some  angles  are  dominant.  The 
calculations  of  reflection  and  refraction  coefficients  as  presented  by 
Muskat  and  Meres  separate  the  problem  into  two  parts  by  considering  the 
incoming  compressional  and  shear  waves  individually.  The  resulting 
expressions  for  reflection  and  transmission  coefficients  are  of  sufficient 
complexity  to  require  computer  calculation.  Such  calculations  have  been 
carried  out  that  were  necessary  to  evaluate  the  Nevada  Test  Site.  They 
show  that  only  one  boundary,  that  between  WL  and  AL,  has  a  significant 
mode  conversion.  In  terms  of  energy,  31  percent  goes  into  the  reflected 
compressional  wave,  64  percent  is  transmitted  through  as  a  compressional 
wave,  and  5  percent  is  a  transmitted  shear  wave  yielding  a  total  of  100 
percent  from  the  incident  compressional  wave. 

3.6  MODEL  OF  TUNNEL  RFSPONSE 

An  estimate  of  the  amplitudes  of  reflected  waves  returning  to  the 
surface  from  the  tunnel  is  required.  Two  feasible  methods  of  obtaining 
the  estimate  are  a  diffraction  model  which  considers  the  tunnel  as  a 
wave  scatterer,  and  a  conical  wave  model  in  which  boundary  conditions  at 
the  tunnel  walls  are  considered  and  the  amplitude  of  emitted  conical 
waves  are  calculated. 

Investigation  of  these  two  approaches  indicated  that  the  conical 
wave  model  was  more  useful.  This  model  yields  values  of  displacement 
for  any  point  in  a  geophone  array. 
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„  _PiscJJ88ion  °f  conical  wave  models  can  be  found  in  Seismic  Waves, 
gadiation  Transmission,  and  Attenuation  by  White.16  These  models, 
developed  primarily  to  investigate  the  wave  propagation  along  cylindrical 
bore  holes,  are  based  on  the  assumption  that  incident  waves  on  the  bore 
hold  have  a  wavelength  long  compared  with  bore  hole  diameter.  The  cylin¬ 
drical  symmetry  of  the  bore  hole  then  determines  the  wave  modes  which 
it  can  emit  when  the  walls  are  displaced.  Both  compressional  and  shear- 
type  waves  are  possible.  Axial  symmetry  permits  describing  the  wave 
field  with  two  potentials. 

The  potential  equations  are 


+ 1  i£  +  i!i  _  i_  i± 

3r2  r  3r  a  2  2  a*2 

°  r  3z  a  3t 


a  2  r  3r  2,2 
3r  r  3  z 


1 

7 


Where  a  is  compressional  conical  wave  speed 
@  is  shear  conical  wave  speed 
r  is  radial  distance  from  center  of  tunnel 
z  is  distance  along  tunnel  axis. 

NOTE:  It  is  assumed  throughout  that  the  conical  wave  speed  does  not 
differ  significantly  from  plane  wave  speed.  Similarily,  spherical  wave 
speed  is  the  seme  es  for  plane  waves. 

The  displacements  can  be  obtained  from  the  potentials  as  follows: 


u 


r 


3$  _ 

3r  3z 


u 
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1£  +  It  +  i 

3z  3r  +  r 
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The  potentials  obtained  by  solving  the  above  equations  for  the  case 
8  <  a  <  | c |  where  c  is  phase  velocity  are 


♦c  -  ^Hq1  0»r)  +  A2Hq2  (mr)j  e~ilz  eio,t 
H>c  -  [ BjHj1  (kr)  +  B^j2  (kr)  j  e  ilz  eiu,t 

9  9  1/2  22 

£  ■  d)/c;  m  ■  to  (1/a  -  1/c  )  ;  k  "  o;  (1/8  -  1/c  ) 


where  u>  -  2n  (frequency) 
c  is  phase  velocity 
H  indicates  a  Hankel  function 

Before  applying  these  equations,  it  is  necessary  to  evaluate  their 
validity  for  the  Nevada  Test  Site.  In  addition,  it  is  useful  to  simplify 
the  relationships  as  much  as  possible.  The  compress ional  wave  speed  in 
the  rock  surrounding  the  tunnel  is  19200  ft/sec  and  the  wavelength  is 
38.4  feet  at  500  Hz.  At  200  Hz,  the  wavelength  is  96  feet.  Tunnel 
diameter  can  vary  from  20  feet  to  30  feet  so  that  at  maximum  frequency, 
the  wavelength  approaches  tunnel  dimensions;  however,  with  the  vibrator 
source,  the  frequency  can  be  adjusted  to  increase  the  ratio  between 
wavelength  and  the  tunnel  diameter  to  any  desired  value. 

Of  particular  interest  are  ratios  of 


uia/a  <  1.0 


where  a  is  tunnel  radius.  For  a  30  foot  tunnel  a  «  15  feet,  and  with 
a  -  19200  ft/sec  this  implies  a  frequency  of  200  Hz  or  less.  Under 
these  conditions,  the  reflected  wave  energy  is  concentrated  in  the  com- 
pressional  wave  for  decreasing  values  of  the  phase  velocity.  At  the 
Nevada  Test  Site,  the  layering  is  such  that  angles 


Y  »  arcsin  (a/c) 
P 


\ 
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will  not  be  less  than  60  degrees  and  the  reflected  congressional  wave 
will  never  contain  less  than  75  percent  of  the  incident  wave  energy 

nnn0^r?r  iC  ls  Posable  to  omit  the  mode  conversion  and 

TMold?  th?  comFressional  wave  contribution  to  the  displacements 

This  simplification  is  not  mandatory,  but  only  a  convenience for  caJcuia- 

con aid-  at  the  NeVada  Teat  Sire  a  more  detaUed 

consideration  of  conical  shear  waves  can  be  conducted. 

The  displacements  reduce  to 


u 


r 


3r 


u 

z 


i! 

3z 


The  calculational  procedure  is 
tunnel  response* 


as  follows  for  determination  of 


or  rh.1,  :5n!ralVU  is  only  necessary  to  consider  a  source  located 
charge  ut  T  for  -  «»U  explosive 
v^e8wlth  pot«"l.f !r8“  8°Ur“  °f  “  Bpherl“l  comPr*ssionai 


8 


.  „  i(kr*-wt) 
s  r* 


where  k  is  the  wave  number  (u/a) 

r  is  radial  distance  from  source* 
Theddisplacement  ufl  Is  obtained  in  the  usual  way 


u 


r 1 


A  w 

a  „  -  S  O 

3rf  r*  a 


SlSV"!’  Ur’  iS  the  knOVm  input  ValUe  and  rorresP°nd8  to  a  fixed 
Hence,  Ag  is  determined 


A 

s 


(u)Q/a) 
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Now  ur«  can  be  determined  for  any  r',  -nd  then  multiplication  by  the 
attenuation  factor  exp(-ar’)  yields  t>-  a  final  value  at  the  tunnel 
boundary. 

For  a  local  region  in  the  vicinity  of  the  tunnel,  it  is  possible 
to  approximate  the  displacement  of  the  spherical  wave  by  a  plane  wave. 
The  potential  for  the  plane  wave  can  be  written. 

-iiz  iwt 

d>  =  A  e  e  e 
P  P 


The  displacements  are  defined  by 


u 

x 


14E 

3x 


3z 


The  relationships  between  u^,  and  the  plane  wave  displacements 


arc 


u  B  u  ,  cos 9 
x  r 


u 

z 


u  ,  sin0 

r 


where  0  is  angle  to  the  vertical. 

The  value  of  Ap  is  determined  as  in  the  spherical  wave  case  and 
the  approximate  plane  wave  is  completely  defined. 

The  incident  plane  wave  upon  arrival  at  the  tunnel  displaces  the 
tunnel  walls.  This  displacement  will  cause  the  formation  of  conical 
waves  and  surface  waves  at  the  tunnel  walls.  At  angles  of  incidence 
near  to  vertical,  coupling  to  conical  waves  will  be  predominant.  This 
is  similar  to  the  reflection  at  a  plane  boundary  where  coupling  to  sur¬ 
face  wave  decreases  with  increasing  angle  of  incidence.  The  establish¬ 
ment  of  completely  symmetrical  conical  waves  is,  of  course,  impossible 
with  an  incident  plane  wave;  however,  an  approximation  is  possible  in 
the  illuminated  zone  and  at  large  wavelength-to-tunnel-diameter  ratios. 

The  two  potentials  $p  and  cf>c  can,  of  course,  only  be  matched  at 
a  single  point  because  the  plane  wave  amplitude  decreases  with  increased 
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SLVrtf  “ave  dl"lnl8hes  ”lth  ‘-casing  radius  from  the 


$ 


P 


A  _iOX  -ill  iut 

Ap  e  e  e 


A1Hq1  (mr)  e"1£z  elut 


♦ 

c 


Equating  real  and  Imaginary  parts 


Ap  cos  ox  -  AjJpOnr) 


Ap  sin  ox  »  AjNpOnr) 

where  r  is  tunnel  radius  and  at  the  tunnel  wall  r  •»  a. 

The  reflection  coefficient  A^  can  be  obtained  as 

a2  Q  ma  H11  (ma)  +  Hq1  (oa) 

A1  Q  ma  Hx2  (oa)  +  HQ2  (oa) 

where 

2 l2  ka  HQ2  (ka) 

(k2  +  l2)  Hx2  (ka)  _ 


.  1  (k2  - 

a2  (k2  -  fc2)2 


Once  A2  has  been  determined, 
is  defined 


the  outgoing  compressional  conical  wave 


(conical  outgoing)  "  A2  H0 


(mr)  e'ltz  s1”* 


For  large  r  (i.e 
Hankel  function  can  be 


,  large  radial  distance  from  the  tunnel) ,  the 
approximated  as  an  exponential. 
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.  .1/2  -l(mr-it/4)  -Hz  eiwt 

.  =  A,  (2/mr)  e 

c (conical  outgoing)  * 


The  partial  derivative  of  this  potential  with  respect  to  r  yields 
the  displacement  due  to  the  outgoing  conical  wave. 


1/2 


1/2  iu/4  e"1^  -Hz  tot 

A,  m  e  T7T  e 


To  evaluate  ur  at  the  boundary  between  AL  and  UAL,  it  must  be  multiplied 
by  the  attenuation  exp(-ur). 

Calculations  of  the  type  described  above  have  been  carried  out  to 
estimate  tunnel  response.  It  Is  apparent  from  the  mode1  that 

response  will  be  obtained  when  the  seismic  source  is  as  close  as  P 
tothetunnel.  At  the  Nevada  Test  Site  this 

sion  to  shear  waves  for  a  geophone  spread  parallel 
3.7  SURFACE  RAYLEIGH  WAVES 

Surface  Rayleigh  waves  are  expected  to  be  grated  whenever  a 
vibrator  Is  operated  at  the  surface.  These  waves  are  usually  of  greater 
amplitude  than  any  other  on  a  seismic  record. 

For  normal  reflection  records,  the  source-geophone  ^ad  i« 

with°the  passage  lY^e  lurf^ilS^w^" rom^a  tranjient^input 
pulse.  With  holographic  viewing,  longer  input  signals  easeful 

sources.  When  layering  occurs  in  the  propagation  medium  such  as  a 
the  Climax  Stock  area,  the  Rayleigh  wave  may  be  dispersive.  1 

characteristic  Rayleigh  dispersion  may  be  found  in  Blast ic_Wa - 

Layered  Media  by  F.wing,  Jardetsky,  and  Fress. 

Generally,  If  the  wavelength  of  the  Rayleigh  wave  Is  shore 
to  the  layer  thickness,  the  layer  appears  as  a  half-space.  The  Rayleigh 

wave  will  not  be  dispersive  and  the  Rayleigh  jjJPt£*foo  Haj 

of  the  shear  velocity.  For  the  high  frequencies  (200  H to 500  Hz) 
to  be  used  in  holography  experiments .  the  shear  “evelength  will  be  very 
short  due  to  the  very  low  shear  ve.ocicy  o  e  ^ 

(860  ft/sec).  The  equivalent  wavelength  at  200  Hz  is 
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thickness  of  the  near-sur f ace  ia/or  is  20  feet  so  that  the  wavelength 
is  only  a  fraction  of  the  thickness.  '.once,  it  is  expected  that  the 
Rayleigh  wave  will  be  nondispvrsi ve  and  have  a  wavelength  of  approximately 


4  feet.  Layering  variations  along  ‘hi 
analysis  and  suitable  cori net  tons  >  .m 
gained  from  the  refraction  surv  y.  N 
culating  the  variations  in  Kty.*  1  ^  «*. 

known. 

Attenuation  N  t 
be  accomplished  by  a .  c  i 
geophone  emplac  me:.t  s.n  •  . 

Test  Site  if  requited. 

3.8  TUNNEL  DES i ki  :> b.  '  E  ‘  .  . 

TRANSMISSION 

A  des tressed  c  aie  1.  I  i  ...  j  ,  >  ' 

drive  tunnel  comples.  nhls  o  •  . ,siJ 
rock  which  has  variant*  ml  -.r 

average  depth  1  oi  the  z  *■:.«.  j  i  , 

of  fracturing  of  the  .  ..  k  ir  /  K  •  * 

The  primal  y  i ^  i  .  '•  i 

the  incident  wave  Iw  .;t 

appear  larger  to  t._ 
occur .  Data  \  i  . . 

degree  v-f  f :  »tt  .r  .„ 
tressed  zone  r.inbc  t 
is  about  11000  ft/  Cv.  Ih 
be  calculated.  An  it  i 

from  an  incidert  c^ir.^r^  : 

CASE:  Incident  cc 

with  vertical. 

Re f  lec to  i  a.\  ...  j  * 
Refle.u  i  .  * 

Transmi  t  L  .1  -  .  . 

Transmitted  sht ai 
Total  In  id  t  t  j  ' 

Therefore  ,  mode  c  >r.  i  cm  c.  ;  «.• 
will  vary  according  to  ii,l.lc«.,u  .  ;.Kfc 
of  compressional  energy  is  r.-jk  I*  j  ; 
The  incident  wave  uill  i.w*  ..l.m 
tunnel  wall. 


geoph^ne  spread  can  modify  this 
made  a;cording  to  information 
J.ffic  ilty  is  anticipated  in  cal- 

•  ,7.v  de  sinc<-  the  theory  is  well 

.i  (\.it;nuou;  source  ran  often 
components  from  a  multi- 
r.  an  will  k:  e  used  ;»•;  the  Nevada 


.  k;  t.:.g 

to  exist  si.i  round  i  ng  the  Pile- 
i  is  a  j  w  velocity  none  of 
.  .  f rem  2  feet  to  20  feet.  An 
II,  low  v.  Kelly  is  the  result 
•  Miing  the  excavation  process. 

.ai  z  no  i  ^  its  (  i  fort  upon 
r.  first,  the  tunnel  can 
r  : ,  mode  conversion  can 
kt  ah...  l  *  i  Cc  *i  1 :  .  on  between 
’  ,  .  ;  .  .  1  u  'he  ces- 

it  a;,  v  rail  average 

v  ...ne  v  a  w.i  e  prupigation  can 
r?.  ion  ic.  L  different  modes 
r  nait.u  be  low. 

Uve  .l  an  ..nale  of  30  degrees 

O'? 

29.4 

7H 

9% 

ioo  x: 

.  .  a  a  ii,  l  r  ib  j t  i  >n  '  f  energy 
...  t is  case,  no  ref  lection 
a  .1  due  to  tne  destressed  zone, 
s  a.  equ^ntly  r^fler  tod  by  the 
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The  Piledriver  data  indicate  thatf fracturing  is  a  major  factor 
inrreducing  wave  velocity.  There  is  therefore  some  concern  as  to  the 
effects  of  the  nuclear  detonations  which  have  occurred  near  the  tunnel 
complex.  Figure  3-7  is  a  cross  section  of  the  Hard  Hat  detonation  site 
and  Figure  3-8  shows  the  extent  of  crushed  and  collapsed  zones.  Based 
upon  geologic  observation,  no  explosion-produced  fractures  are  reported 
from  the  15U0  entrance  shaft  and  extending  280  feet  toward  the  U15a 
shot  point.  The  actual  effect  upon  elastic  wave  propagation  can  only 
be  inferred,  since  no  post-shot  geophysical  investigations  have  been 
carried  out. 

At  the  Piledriver  site,  the  tunnels  vere  resurveyed  after  the 
detonation  and  permanent  displacements  were  noted  at  all  points. 

Displacement  at  1000  feet  is  0.86  foot 

Displacement  at  1500  feet  is  0.38  foot 

These  measurements  were,  of  course,  made  on  the  destressed  zone  and 
the  validity  when  extended  into  the  higher  velocity  rock  surrounding 
the  tunnel  is  uncertain.  Of  greater  importance  and  significance  are 
bore  hole  data  from  the  U-15.01  PS  IV  drill  hole  Survey  of  Piledriver 
Results ,  by  I.  Y.  Borg.20  This  drill  hole  penetrated  into  the  cavity 
formed  by  the  Piledriver  event.  Permeability  data,  which  are  related 
to  degree  of  rock  fracturing,  indicate  a  zone  of  increased  permeability 
with  a  maximum  vertical  extent  of  1032  +  36  feet  above  shot  level  or 
approximately  500  feet  from  the  surface.  The  Climax  Stock  has  many 
fracture  planes  and  these  have  mostly  been  healed  with  secondary  minerals. 
The  Piledriver  detonation  has  most  certainly  opened  large  numbers  of 
these  fractures  at  least  to  a  distance  of  1050  feet  from  the  shot  point. 
The  effect  beyond  1050  feet  is  uncertain,  but  it  must  in  any  case  be 
assumed  that  some  velocity  and  attenuation  changes  have  taken  place. 

It  is  also  quite  likely  that  these  opened  fractures  are  now  filled  with 
perched  water.  For  this  reason  and  for  layering  corrections,  the  pre¬ 
liminary  refraction  and  reflection  tests  are  recommended. 
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Figure  3-7  -  Generalized  Cross  Sectio.:  of  1500  Tunnel 
Showing  Relation  of  Tunnel  to  Shot  Point 
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Figure  3-8  -  Generalized  Hap  of  1500  Tunnel  Showing  Crushed 
and  Collapsed  Zones  Due  to  the  U  15a  Explosion 
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SECTION  4 


DESIGN  OF  UNDERGROUND  VIEWING  SYSTEM  - 
SEISMIC  HOLOGRAPHY  EXPERIMENT 


In  this  section  the  design  of  the  underground  viewing  system  -  or 
more  accurately  -  the  design  of  the  seismic  holography  experiment  is  dis¬ 
cussed,  Based  on  the  analysis  of  the  selected  field  test  site  geology 
described  in  Section  3,  calculations  are  made  to  predict  signal  strengths 
present  at  the  surface.  These  signals  include  those  from  the  targets ,  as 
well  as  reflections  and  refractions  from  reflecting  interfaces.  Once  the 
signal  strengths  are  known,  the  following  information  can  be  obtained: 
the  maximum  frequency  (from  a  vibrator  energy  source)  that  can  be  used  to 
ensonify  the  targets,  the  maximum  spatial  frequency  (minimum  fringe  spac¬ 
ing)  across  the  aperture,  the  number  of  sampling  points  required  in  the 
holographic  aperture,  and  the  size  and  position  of  the  aperture.  The 
latter  can  be  determined  by  considering  the  target  signal  strengths  and 
the  topographical  features  of  the  site.  Data  acquisition  with  different 
types  of  seismic  excitations  (CW  and  pulsed  CW  using  a  vibrator,  and  im¬ 
pulses  using  explosives)  are  considered.  Finally  the  equipment  and 
instrumentation  requirements  are  given. 

4.1  CALCULATED  AMPLITUDES  OF  RECEIVED  SIGNALS  BASED  UPON 

THE  COMPOS ITi  MODEL  OF  THE  NEVADA  TEST  SITE 

The  vibrator  source  is  described  and  the  expected  ground  displace¬ 
ments  are  given  in  Section  5.  With  these  values  of  displacement  as  in¬ 
puts,  expected  amplitudes  from  important  wave  paths  have  been  calculated. 
The  values  obtained  consider  the  tunnel  response,  geometrical  spreading, 
attenuation,  energy  loss  through  mode  conversion  at  layer  boundaries, 
the  destressed  zone  surrounding  the  tunnels,  and  coupling  of  the  incident 
wave  into  an  outgoing  conical  wave. 

A  complete  calculation  is  carried  out  as  follows:  At  500  Hz  the 
surface  displacement  at  the  source  is  5  x  10”^  in.  The  wave  enters  the 
ground  almost  vertically  downward  at  an  angle  of  incidence  from  normal 
of  30  degrees  for  the  Pil^driver  case  (see  Figure  4-1).  The  wave  prop¬ 
agates  downward  as  a  spherical  wave  and  attenuates  due  to  geometrical 
spreading  as  1/R,  In  addition,  the  wave  is  attenuated  at  a  rate  which 
is  determined  by  its  frequency.  The  attenuation  factor  is  exp(-aR), 
where  a  corresponds  to  attenuation  factors  previously  calculated.  These 
factors  combine  to  determine  the  wave  amplitude  in  the  near-surface  layer 
just  before  crossing  into  the  weathered  layer. 

II  ^r  (incident)  -aR 

R(final)  R  e 
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SOURCE  GEOPHONE  DISTANCE  1,380  FT 
OURCE  NEAR  SURFACE  LAYER  20  FT.  GEOPHONE 

- X - ^.WEATHER^FT^Tiil^lT"""""— 

\  *1 TEREO  LAYER  JOoct  l. 


UNAl  T&REO  lave  r  1 . 1 00  FT. 


TU^PlCt  DlAMtten  TDFt. 


SCALE:l 


!  200  FT. 


Figure  4-1  -  Herd  Her  Slt.,  Ray  Paths  Sample 

r(incident)  "  x  10  inch 
R  (thickness  of  layer)  ■  20  feet 

_  SM  D 


e~aR  -  0.705 


^r(final)  *  1*76  x  10  ®  inch 


Hi  ^  °5  partlculsv  interest  is  the  vertical 

t(h^MQTSlnue  thl8  18  the  c°“P°nent  measured aJvcomponent  of  displacement 
L  the  incident  angle  is  3  degrees  and  U  8eoPhones.  For 

To  determine  the  effect  nf  n  k  X  r* 

mission  coefficients  are  calculated  the  reflection  and  trans- 

“ 18  2* 


4-2 


In  this  case,  u  “  2nf  *  31A0  and  Vr  ■  2150  ft/sec. 


I  =  1.21  x  10"6 

I2  -  ^  <V  ■  1.51  *  10'6 

Vertical  displacement  in  the  weathered  layer  becomes 

?  -  1.32  x  10'6  inch 


U 


2  V 


The  angle  of  incidence  for  the  weathered  layer  is  5  degrees,  so  that  Ux 
is  approximately  equal  to  Ur. 

The  wave  amplitude  is  again  corrected  for  geometrical  spreading 
and  attenuation  factor  a  corresponding  to  the  weathered  layer. 

U  ...  ,,  -  1.16  x  10~7  inch 

r (final) 

A  new  value  of  the  incident  potential  amplitude  is  calculated 

»-7 


Xl“ 


U  V 
x  t 


1.33  x  10 


where  V  corresponds  to  congressional  wave  speed  in  the  weathered  layer. 

The  potential  amplitude  ratio  for  the  altered  layer  is 
I2/!i  -  1.56.  Thus, 


I2  -  2.07  x  10 


-7 


U  =•  A. 58  x  10 
x 


-8 


Correcting  Ux  for  geometrical  spreading  and  attenuation  yields 


U  ,v  -  8.62  x  10 

x(final) 


-8 


The  new  potential  is 

I  =  3.53  x  10 


A-3 


For  the  boundary  between  altered  layer  and  unaltered  half-space 


—  -  1  '’4 

■i 


I2  -  4.37  x  10 


-8 


U  -  7.08  x  10 
x 


-9 


Correcting  for  geometrical  spreading  and  attenuation 

U  ...  -  1.26  x  10~10 

x  (final) 

Therefore,  the  potential  wave  amplitude  Is 


1^  -  / .78  x  10 


-10 


In  order  to  calculate  the  tunnel  response  It  Is  necessary  to  deter¬ 
mine  the  potential  of  the  outgoing  conical  wave.  First,  the  congressional 
plane  wave  potential  Is  matched  to  the  potential  of  a  conical  wave  at  the 
boundary  of  the  tunnel.  This  matching  can  only  occur  on  the  ensonifled 
side  of  the  tunnel  and  along  the  center  line  of  the  tunnel.  The  contri¬ 
bution  of  Ux  to  a  radial  displacement  of  the  tunnel  wall  drops  off  rapidly 
along  either  side  of  the  center  line.  The  emitted  conical  wave  will 
therefore  not  be  symmetrical. 

Matching  Is  carried  out  for  real  and  Imaginary  parts 

Real  part:  cos  (mx)  -  JQ  (ma) 

m  -  -  1.62  x  10"1 

P 

Ax  -  -1.05  x  10~9 

Imaginary  part:  ^  sin  Ai  N0  (ma) 


A^  -  1.46  x  10 


Using  ihe  previously  developed  expressions  for  A  and  A2/A^, 
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Q  »  -0.318  -  1  (0.337) 


A2/Al  “  -°*033  +  1  (1,65) 


The  real  or  imaginary  part  of 
conical  wave  is  positive  and  real, 
this  condition. 

iut 

e  -*■ 

ilz 

e 

imr 

e 

in/4 

e 


must  be  chosen  so  that  Ur  of  the 
The  following  are  chosen  to  satisfy 

i  sin  ut 
cos  Hz 

cos  mr 
►  cos  n/4 


The  amplitude  is  determined  when  sin  u>t  ■  1. 

cos  £z  ■  0.809 


cos  mr  ■  -0.309 


Hence 


1.46  x  10~9 


is  chosen  to  yield  positive. 


U  ■  4 . 08  x  10  inch 
r 


Including  attenuation 


^r(final) 


;.07  x  10 


-13 


inch 


Next,  the  potential  is  calculated  as  if  the  emergent  conical  wave  were 
an  approximation  to  a  plane  wave. 


-  3.13  x  10 


-12 


For  the  boundary  of  unaltered  layer  into  the  altered  layer,  "  0*748 


I2  -  2.34  x  10~12 

U  -  5.72  x  10~13 
x 
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Including  conical  wave  spreading  and  attenuation. 


Ux (final) 


-  2.96  x  10"13 


The  potential  becomes 


1.21  x  10 


-12 


For  the  boundary  of  altered  layer  and  weathered  layer  I2/Ii  -  0.337. 
tiSnaJieldshe  CalCUlatlon  as  above  lncluding  wave  spreading  and  attenua- 


U 


x(final) 


-  1.51  x  10 


-13 


1^  ■  1.73  x  10 


-13 


I2^I1  *or  tbe  weathered  layer,  near-surface  layer  boundary  is  0.742 

,-13 


I2  •  1.28  x  10 


Wave  spreading  and  attenuation  yield  a  surface  value  of 


U 


x(flnal) 


1.3  x  l(f13  inch 


a<  <;°n8ldera'lon  has  so  far  been  given  to  the  destressed  zone  sur- 
th«  tunnel-  The  zone  will  cause  mode  conversion  and  calculation 
°f  indicate  that  only  71.3  percent  of  the  incident  wave  energy  will 

reech  the  tunnel  through  the  destressed  zone.  Upon  re-emerging,  the  wave 
will  lose  12.3  percent  of  wave  energy  to  a  trapped  mode  surfouSding  the 

10M  yleld‘ *  r“ul“"t  ■“*-  lil 

T.  m  *  compilation  of  expected  signal  returns  is  shown  in  Table  4-1. 
Theue  include  the  response  of  Hard  Hat  (Figure  4-2)  as  well  as  Piledriver 
tunnels  to  various  frequency  inputs, 

KXt  86en  fc!?at  th68e  re8Ponsea  are  quite  low  and  consideration 

dPtirM  8lVfn^J°  the  p088iblllty  that  other  8ignals  may  interfere  with 
detection  of  th  “.unnel  signal. 

Will  ,-pn!!e!:ent  8i8naltreturns  from  layering  and  surface  wave  interference 
Wiil  tend  to  use  up  the  dynamic  range  of  the  recording  system,  particu- 

siblP  the  J"put,ls  a  continuous  wave  (CW) .  An  investigation  of  pos¬ 

sible  strong  reflections  and  refractions  has  been  undertaken.  The  ? 
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Table  4-1  -  Wave  Amplitudes  for  Piledriver  and 
Hard  Hat  Test  Sites 


Input  displacement  from 
vibrator 

Tunnel  Response,  Hard  Hat 
Source-geophone  distance 
695  feet 

Tunnel  Response,  Piledriver 
Source-geophone  distance 
1380  feet 

P  refraction  at  695  feet 

P  refraction  of  1380  feet 

Rayleigh  Wave  at  695  feet 

Rayleigh  Wave  at  1380  feet 

P  reflection  (15  degrees 
from  vertical) 


500 

Hz 

200 

Hz 

100 

Hz 

5  3t 

10~5 

2.5  x 

10-3 

1  X 

io”2 

5.3  x 

io"12 

4.7  x 

10“9 

5.4  x 

ic-8 

9.7  x 

io-14 

8.2  x 

IQ"10 

1.6  x 

10-8 

7.7  x 

10-11 

1.2  x 

10”8 

1.1  X 

io-7 

7.7  x 

10"12 

4.4  x 

10~9 

8.6  x 

io-8 

7.1  x 

10-8 

3.6  x 

10~6 

1.4  x 

io-5 

3.6  x 

10-8 

1.8  x 

io“6 

7  x 

io”6 

1.6  x 

10‘8 

8.1  x 

10-7 

3.2  x 

io-6 

strongest  reflections  from  the  reflecting  layer  will  occur  near  the 
source.  P-waves  entering  the  ground  at  15  degrees  from  the  vertical  and 
reflected  from  the  WL-AL  interface  will  give  a  ground  displacement  ampli¬ 
tude  as  tabulated  below. 

Reflection  (15  degrees  from  vertical,  45  feet  from  source) 

500  Hz  -  1.6  x  IO-8  inch 
200  Hz  =  8.1  x  10  ^  inch 


100  Hz  ■  3.2  x  10  8  inch 


P-waves  entering  the  ground  at  an  angle  greater  than  15  degrees 
(the  critical  angle)  will  be  refracted  rather  than  reflected.  The  strong¬ 
est  returning  body  wave  due  to  this  refraction  is  the  PiP2p3p2pl  refrac¬ 
tion.  This  wave  enters  as  a  compressional  wave,  is  refracted  by  the  WL-AL 


SOURCE  GEOPHONE  DISTANCE  690  FT. 


Figure  4-2  -  Piledriver  Site,  Ray  Paths  for  Sample  Calculation 


boundary,  travels  with  compressional  wave  speed  in  AL,  and  then  returns 
to  the  surface  as  a  compressional  wave.  The  discovery  that  this  is  the 
strongest  body  wave  return  is  not  surprising,  since  it  is  the  prominent 
wave  mode  reported  in  the  V.'la  Uniform  reports.13  Amplitudes  for  this 
refraction  are  tabulated  below. 


Refraction  at  695  feet 

500  Hz  =  7.7  x  10"11  inch 

200  Hz  =  1.2  x  10"8  inch 

100  Hz  =  1.1  x  10"7  inch 

Refraction  ct  1J80  feet 

500  Hz  r  7.7  x  10"12  inch 

200  Hz  =  4.4  x  10  8  inch 

100  Hz  =  8.6  x  10'8  inch 
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These  amplitudes  are  not  sufficiently  great  to  exceed  the  dynamic  range 
of  the  recording  system  and  therefore  will  not  mask  the  tunnel  response. 

The  surface  wave  at  the  higher  frequencies>  greater  than  100  Hz, 
is  essentially  not  dispersive  and  can  be  attenuated  by  a  geophone  group. 
Methods  for  estimating  the  effective  attenuation  can  be  found  in  A  Graph¬ 
ical  Method  for  Computing  Geophone  Group  Response t  by  Verma  and  Roy . 2 1 
For  a  continuous  sinusoidal  input,  the  normalized  response  Rn  of  n  single 
geophones  of  equal  sensitivity  in  a  linear,  equally  spaced  array  is 


sin  (n  tt  D/A) 
n  n  sin  (tt  D/A) 


where  A  is  the  wavelength  and  D  is  the  distance  between  consecutive 
geophones . 


For  geophones  of  equal  sensitivity,  the  effective  attenuation  is 
strongly  dependent  upon  the  exact  placement  of  the  geophones.  This  crit¬ 
icality  can  be  reduced  by  adjusting  the  relative  sensitivity  of  the  geo¬ 
phones  in  the  group.  As  an  example,  a  six-element  array  with  spacing 
varying  from  A/1.5  to  A/4  will  give  a  normalized  response  of  0.06  or  less 
to  a  wave  with  wavelength  A. 


Determination  of  an  optimal  spacing  and  sensitivity  can  be  estab¬ 
lished  by  summing  terms  of  the  form  Fn.  Estimation  of  Rayleigh  wave  amp¬ 
litudes  is  difficult  for  near-surface  materials.  However,  for  low  fre¬ 
quencies,  Jolly  and  Mifsud22  have  reported  that  amplitude  varies  roughly 
as  the  reciprocal  of  shot  detector  distance  for  the  vertical  component. 
This  decay  in  amplitude  is  the  combined  result  of  wave  spreading  and  at¬ 
tenuation.  At  frequencies  above  100  Hz,  the  amplitude  is  expected  to 
decrease  more  rapidly  due  to  higher  attentuation.  A  suitable  attenuation 
factor  has  not  been  determined  for  surface  waves;  therefore,  amplitude 
will  be  calculated  using  the  inverse  distance  relationship  only.  It  is 
recognized  that  this  is  only  an  upper  limit  to  the  expected  amplitude. 

Rayleigh  wave  amplitudes  at  700  feet 
500  Hz;  7.1  x  10  ®  inch 
200  Hz;  3.6  x  10  ^  inch 


100  Hz;  1.4  x  10  inch 
Rayleigh  wave  amplitudes  at  1400  feet 

.-8 


500 

Hz; 

3.6 

X 

10 

inch 

200 

Hz ; 

1.8 

X 

IQ'6 

inch 

100 

Hz ; 

7.0 

X 

io-6 

inch 
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As  pointed  out  previously,  the  use  of  multiple  geophones  can  reduce 
the  effective  amplitude  of  Rayleigh  waves  by  a  factor  of  0.06. 

Effective  Rayleigh  wave  amplitudes  at  700  feet 


500 

Hz; 

4.3 

X 

io-9 

inch 

200 

Hz; 

2.2 

X 

1 

O 

•H 

inch 

100 

Hz; 

8.4 

X 

1 

O 

•H 

inch 

Effective  Rayleigh  wave  amplitudes  at  1400  feet 
500  Hz;  2.2  x  10-8  inch 
200  Hz;  1.1  x  10"7  inch 
100  Hz;  4,2  x  10  7  inch 

The  Rayleigh  wave  amplitudes  will  be  greater  than  those  listed 
above  nearer  the  source.  The  maximum  value,  occurring  a  few  wavelengths 
(for  example  35  feet)  from  the  source,  will  be  approximately  20  times 
greater  than  the  amplitudes  listed  above.  This  indicates  that  the  dynamic 
range  requirements  may  be  exceeded  at  geophones  near  the  source,  but  this 
should  not  be  a  serious  problem  since  the  a  "hole"  in  the  aperture  will 
not  drastically  affect  the  reconstruction. 

fi  Table  summarizes  the  displacement  amplitudes  of  the  various  re¬ 
flections,  refractions,  and  Rayleigh  waves  calculated  in  this  section. 
According  to  this  table,  the  strongest  signal  will  be  due  to  the  Rayleigh 

th^de  the  raCJ°  °f  the  RayleiSh  wave  signal  to  the  signal  from 

the  deep  Piledriver  tunnel  will  be  within  the  dynamic  range  of  the  record¬ 
ing  equipment  at  geophones  greater  than  200  feet  from  the  source.  If  the 
Rayleigh  wave  is  attenuated  by  using  multiple  geophones,  the  ratio  is  re¬ 
duced  by  greater  than  a  factor  of  10. 

4.2  HOLOGRAPHIC  APERTURE  DESIGN 

The  signal  strength  calculations  of  Section  4.1  indicate  that  the 
maximum  excitation  frequency  that  can  be  used  to  ensonify  the  deep  (1400 

from\-hUnnh1nS  2  “Z;  F“rther>  extrapolating  between  the  signal  strength 

hr  Sh®JloW  (78°  feet)  tunnel,  frequencies  as  high  as  300  Hz  might  be 
used.  For  the  purpose  of  the  report,  however,  a  vibrator  excitation  fre¬ 
quency  of  200  Hz  will  be  assumed.  If  initial  field  reflection  results  in¬ 
dicate  that  either  higher  or  lower  frequencies  can  or  must  be  used,  the 

accordingly^  hol°8raphic  aPerture  (geophone  spacing)  can  be  adjusted 

SinrP  Il!nlUSf  °f  eXploslva®  as  a  seismic  energy  source  is  also  recommended. 
v  h«h  !P  iV6S  are  Cfpable  of  listing  much  more  seismic  energy  than  a 
vibrator  source,  even  in  a  narrow  frequency  band,  quasi-monofrequencies  of 
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greater  than  200  Hz  may  be  obtained.  Whether  frequencies  with  high  inten¬ 
sities  greater  than  200  Hz  can  be  realized  will  be  determined  by  the  fre¬ 
quency  analysis  of  records  in  the  field. 

The  maximum  holographic  aperture  that  can  be  used  is  estimated 
to  be  1500  feet  square.  This  estimate  is  based  primarily  on  the  long 
path  signal  strengths  from  targets  as  calculated  in  the  previous  sec¬ 
tion.  Figure  4-3  illustrates  the  outline  of  the  aperture  recommended 
for  the  experiment.  The  size  of  the  aperture,  and  its  position,  was  also 
influenced  by  the  topography  of  the  area.  Roads  run  parallel  to  the  west 
and  south  sides  of  the  aperture  and  the  topography  becomes  more  rugged 
outside  the  aperture  outline  shown. 

The  seismic  energy  source  (vibrator  or  explosives)  is  shown  in  the 
southeast  corner  of  the  aperture  in  Figure  4-3.  This  position,  near  the 
intersection  of  the  two  tunnels,  was  selected  so  that  both  tunnels  could 


Mir  ooo 


fh/a  ooo 


Figure  4-3  -  Map  Showing  Location  of  Holographic  Aperture 
with  Respect  to  Nevada  Test  Site 
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frin^rifled  C°f ther  more  effectively  to  provide  a  maximum  number  of 

chim5pv«aK  ue  aperture'  and  because  it  is  further  removed  from  the 
chimneys  above  the  two  nuclear  shot  points. 

Mono  n|™ber  oF  fringes  obtainable  across  the  aperture  due  to  reflec- 
wee  J  Zh  ""  dct'mlncd  ^  calculating  the  maximum  dlffer- 

.Ea  .  .E 8Eh  from  a  polnt  on  the  tunnel  vertically  to  the  surface 
and  to  the  edge  of  the  array  (see  Figure  4-4).  Velocity  changes  in  the 

frequenev^f  20oT  T  acc°unted  for*  Fr°“  these  calculations,  at  a 
fle?Mnno  f  20  the  number  of  fringes  across  the  aperture  due  to  re¬ 

tunnel  wilirRivehse^ePftl,nnel  18 T!llghtly  more  than  four‘  The  shallow 
tunnel  will  give  seven  fringes.  The  use  of  higher  or  lower  frequencies 

will  increase  or  decrease,  respectively,  the  nSmber  of  fringes 

na  The  minimum  fringe  spacing  across  the  aperture  will  determine  the 

rX  sepX™  1  8ha°P5T  TCJnf)  case L. 

P^om  Fiwra  4?4  t  ?  "Cd<by  rcfl«tions  from  the  shallow  tunnel. 

owncy  of  200  Ha  E.  T  BPaa1"8  ls  60  fcct  f»r  “  excitation  fre- 
quency  ot  200  Hz.  Since  a  function  is  adequately  sanmled  if  «f  ioa„. 

samples  per  cycle  (fringe)  are  taken.”  a  „axl»u„  gMJho„e  spacl„1ga“ft  Cw° 
spacing  5  ItTe'X  5£SS.~«““  ‘  ^Ue 


MINIMUM 
FRINGE  SPACING 


TUNNEL 

Figure  4-4  -  Determination  of  Number  of  Fringes 
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The  angular  resolution  a  of  a  square  aperture  A  with  an  Illuminating 
wavelength  X  is24 

a  ■  sin  *  (A/A) 

For  this  case,  with  f  *  200  Hz, 

a  »  sin”1  (100/1500)  ■  0.067  radians 

Since  the  two  target  tunnels  are  780  and  1400  feet  deep,  the  minimum  dis¬ 
tance  A  resolved  will  be 


A,/rtrt  _  ,  ■  1400  a  ■  94  feet 

1400  Tunnel 


700  Tunnel 


760 


50  feet 


The  resolvable  distances  are  larger  than  the  tunnel  diameters  (20  feet), 
hence  images  of  the  tunnels  in  the  reconstructions  will  appear  wider 
than  they  actually  are. 


4.3  HOLOGRAPHIC  DATA  ACQUISITION 

Holographic  data  can  be  gathered  using  four  different  types  of 
excitation:  monofrequency  CW,  pulsed  monofrequency  CW,  and  frequency 
modulated  or  "chirp"  pulses  with  a  vibrator  source,  and  with  impulsive 
sources  using  explosives. 


4.3.1  Monofrequency  CW 

Ensonification  of  a  target  area  with  a  continuous  mono- 
Crequency  wave  long  enough  to  establish  steady  state  conditions  is  the 
method  that  usually  comes  to  mind  when  holography  is  considered.  Using 
this  method,  every  sampling  point  in  the  aperture  simultaneously  re¬ 
ceives  sigials  from  every  reflecting  and  scattering  point  and  from  every 
noise  sourte  in  tV  earth  within  its  range  after  steady  state  conditions 
are  establr.shed.  Theoretically,  one  could  instantaneously  sample  each 
point  in  the  array  and  obtain  a  hologram. 

The  disadvantage  of  this  method  in  the  seismic  case  is  that 
strong  reflections  and  refractions  and  signals  from  various  noise  sources 
(wind,  traff.'c  noise,  etc.)  are  present  along  with  weak  signals  from  the 
desired  targets.  Random  noise  can  be  virtually  eliminated  from  the  re¬ 
ceived  signal  using  phase  detection  and  integration  schemes,  but  the 
relative  amplitudes  of  coherent  reflecting,  refraction,  and  scattering 
noises  to  the  desired  target  signal  remain  high.  The  ratio  of  these 
strong  undesired  signals  to  desired  target  signals  approaches  80  dB  near 
the  source  (see  Table  4-1) ,  and  in  the  present  case  is  barely  within  the 
dynamic  range  of  digital  recording  equipment. 
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The  large  amplitude  ratios  of  undesired  strong-to-weak  tar¬ 
get  signals  indicate  that  weak  signal  enhancement  techniques  (described 
in  Section  6)  might  be  used  in  order  to  image  the  target  tunnels  when  a 
continuous  monofrequency  wave  is  used  for  seismic  excitation. 

An  alternate  method  of  defining  the  weak  target  signals  in 
the  presence  of  large  interfering  signals  Is  to  subtract  the  large  signals 
from  the  total  signal.  The  amplitude  and  phase  of  the  Rayleigh  wave  at 
each  geophone  can  be  determined  from  a  pulsed  CW  experiment  (see  Section 
4.3.2).  With  this  information,  the  Rayleigh  wave  can  be  effectively  re¬ 
moved  from  the  total  signal  during  computer  processing.  Strong  signals 
from  the  WL-AL  reflecting  interface  may  also  be  removed  from  the  total 
signal  in  a  like  manner.  The  phase  and  amplitude  of  the  interfering  sig¬ 
nals  must  be  measured  accurately  to  accomplish  this,  however. 

Figure  4-5  shows  the  method  in  which  data  acquisition  is 
accomplished.  The  geophone  signal  is  amplified  by  a  seismic  amplifier, 
mixed  with  a  reference  signal  and  phase  detected,  integrated  to  eliminate 
random  noise,  and  recorded  on  a  digital  recorder  after  analog-to-digital 
conversion.  The  "raw"  received  signal  would  also  be  recorded. 

4.3.2  Pulsed  Monofrequency  CW 

Ensonlf ication  with  pulsed  monofrequency  waves  differs  from 
the  CW  case  in  that  the  pulse  length  is  not  long  enough  to  establish 
stea  !y  state  conditions.  Hence,  signals  received  at  the  sampling  points 
in  the  aperture  vary  with  time,  and  individual  signals  originating  from 
various  reflectors  and  scatterers  and  the  Rayleigh  wave  signal  can  be 
separated  on  a  travel-time  basis.  Random  noise  can  be  reduced  from  the 
received  signals  by  recording  and  "stacking"  (averaging)  the  received 
signal  of  several  pulses.  The  stacking  process  increases  the  signal-to- 
random-nolse  ratio  by  a  factor  of  yN,  where  N  is  the  number  of  records 
averaged. 

The  advantage  of  this  technique  is  that,  although  the  rela¬ 
tive  amplitudes  of  strong  and  weak  signals  are  not  changed  from  the  CW 
case,  many  of  the  undeslred  strong  signals  can  be  discriminated  against 
in  the  time  domain. 


Figure  4-5  -  Holographic  Data  Acquisition 
of  Continuous  Monofrequency  Waves 
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The  method  of  data  acquisition  for  this  case  is  shown  in 
Figure  4-6.  The  geophone  signal  is  amplified ,  converted  from  an  analog 
to  a  digital  signal,  stacked  (time  averaged),  and  recorded  on  a  digital 
recorder. 


4.3.3  Impulsive  (Explosive)  Sources 


The  use  of  an  explosive  to  generate  seismic  waves  has  two 
advantages.  First,  explosives  are  capable  of  radiating  much  more  energy 
than  that  achieved  with  vibrator  sources,  even  in  a  narrow  band  of  fre¬ 
quencies.  Second,  if  the  explosive  charge  is  detonated  below  the  major 
reflecting  layer  (100  feet  deep),  the  reflected  and  refracted  waves  no 
longer  exist  in  the  aperture  plane.  A  large  direct  wave  from  the  source 
is  present.  However,  since  this  wave  arrives  at  the  geophones  prior  to 
reflections  from  the  targets,  it  can  be  discriminated  against  by  time 
gating  or  time  discrimination. 


Because  explooives  radiate  a  wide  band  of  frequencies,  the 
received  signal  must  be  filtered  to  obtain  a  quasi-mono frequency  signal. 
This  is  most  conveniently  done  in  a  digital  computer.  The  minimum  use¬ 
ful  bandwidth  Af  that  can  be  analytically  filtered  depends  on  the  dura¬ 
tion  T  of  the  signal  and  is25 


In  this  case  the  minimum  signal  length  is  expected  to  be  0.2  seconds  re¬ 
sulting  in  Af  >*  0.4  Hz.  For  a  center  frequency  of  200  Hz,  the  bandwidth- 
to-frequency  ratio  is  then  1  to  500. 

The  temporal  coherence  L  of  this  monofrequency  wave  is26 


L 


i  2/f2 

1  c  /f 
4  cAf 
,2 


c 

4Af 


Figure  4-6  -  Holographic  Data  Acquisition 
of  Pulsed  Monofrequency  Waves 
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Z&ZumSA1  2?M°°,  fe‘CtSr  in  era"lte-  thlB  e1-  a  coherence 
jiiLvo.  A  ’  feet*  Thls  length  is  mum  greater  than  any  target  to  be 

mlnation  la  assured?'"  th8"  ^  ale"al  Path  la"et,,i  th“8  Coharart  1U“' 

as  need  for  pulsed"™  tT)  lllTTt  £  77  “**  7  ‘"V8"8 

eral  records  is  not  used.  tlme  avera8ln8  of  sev- 


4. 3. A 


Frequency  Modulated  (Chirp  Pulses) 


lem  is  still  flirlv^11  ?  1  u  ofuhol°8raPhlc  methods  to  the  seismic  prob- 
incornoraJ  *  '  We  therefore  feel  that  it  is  appropriate  to 

J" J®  I?016  ^on^entlonal  ^Pe  of  measurement  into  the  test  program 
methods  tHe  h0l°*raPhlc  aPP— h  to  be  related  to  existing 

j  j  /?ne  relatlvely  common  seismic  approach  is  to  use  a  "''him" 

daremiifrL'a;6-  7  f  a”d  8’"»l<”'  delation  tenures  re'* 

.  t  L1"'  °f  a  «*«“  “««“<>  «  scattered  exponent 

(either"!  ll^S?***  “d>  in  pr3Ctlce’  30  *  geophones  is  Ld 

ceivers)"  ”  “,?"ay  8ynthesized  by  moving  one  or  more  re- 

!r^eenfahc0l°g!aPhy  “*  **  ^hc!nei!!e?Ui^eama85nr°!fff!re1^eriUes 

whereas,  the  te^ral°L2EiI!!^ 

Since  a  source  which  can  be  driven  in  a  ,,ehirnH  mannor  a 

4,3,5  Detection  of  Chimney  Formed  by  Nuclear  Exp1r»Hnn 

u  *he  h°lo8raPhic  experiments  described  in  the  previous  sub¬ 
sections  have  been  designed  primarily  to  detect  the  tunnel  tablets  To 
ensonify  the  chimney  area  formed  by  the  Piledriver  nuclear  exolosion  i-h« 
seismic  aourca  should  ha  moved  to  a  point  In  the  ape^i«  ov«  tha  "ijm- 
y.  It  is  recommended  that  an  explosive  charge  be  used  as  the  enerev 
source.  This  recommendation  It  made  primarily  1„  the  Interest  of  ecoLmy. 

one  line  of  geoX-Y  win  S"10”  aatll“- 

aperture.  Withthis  n„»  ?  j  111  be  used  at  a  time  to  sample  the 

sssssSsswsssr 

ture  have  been  interrogated  During  thi  samPlin8  points  in  the  aper- 
that  the  vibrator  source  be* left  in^he  Lm!"!6  °?®ration*  ic  is  desirable 
coupling  remains  as  constant  as  possible.  Sin^ftt^vJSr:^;"!^^ 


4-16 


will  be  fixed,  the  use  of  the  vibrator  to  ensonify  the  chimney  would  re¬ 
quire  "sweeping"  out  the  entire  array  another  time.  an  explosive 

source  is  used  to  ensonify  the  tunnel,  hovrever,  the  data  can  be  acquired 
during  the  original  aperture  sampling  operations. 

The  frequency  component  of  the  explosive  source  used  to  en— 
sonify  the  cnimney  area  will  be  the  highest  useful  quasi-monofrequency 
determined  from  preliminary  experiments  performed  during  the  field  test 
(see  Section  4.4). 

4.4  AUXILIARY  EXPERIMENTS 

The  calculations  previously  carried  out  to  'jetermine  the  important 
parameters  of  the  Nevada  Test  Site  will  be  supplemented  by  a  series  of 
preliminary  tests.  These  tests  will  be  separated  into  two  operations  as 
described  below.  The  purpose  of  each  test  is  listed  and  interpretation 
of  the  records  will  be  done  in  the  field.  In  this  way,  an  optimum  holo¬ 
graphic  recording  array  will  be  specified.  Preliminary  tests  will  be  run 
with  the  geophone  spread  parallel  to  tunnel  and  directly  over  it.  Both 
Hard  Hat  and  Piledriver  tunnels  will  be  investigated. 

4.4.1  Refraction  Survey 

The  reflection  survey  using  the  vibrator  source  will  be 
made  at  200  Hz  and  500  Hz.  The  vibrator  will  be  moved  to  both  ends  of 
the  geophone  line.  Paper  records  will  be  interpreted  on  site  to  obtain 
local  wave  velocities,  evaluate  local  variations  in  layering,  and  deter¬ 
mine  surface  wave  characteristics.  In  addition,  changes  which  have  oc¬ 
curred  due  to  the  nuclear  detonations  will  be  evaluated. 

4.4.2  Reflection  Experiments 

Using  information  obtained  in  the  refraction  survey,  geo¬ 
phones  will  be  spaced  as  they  will  be  in  the  holography  array.  A  reflec¬ 
tion  survey  will  be  run  using  explosives  emplaced  in  the  altered  layer. 
Shooting  will  be  done  from  both  ends  of  the  line  and  paper  records  ob¬ 
tained.  The  purpose  is  to  determine  the  highest  frequency  which  can  be 
obtained  from  explosives  and  provide  a  strong  signal  to  ensonify  the 
tunnel.  The  explosive  charge  weight  will  be  adjusted  to  enhance  higher 
frequencies.  In  addition,  good  control  will  be  obtained  on  vertical 
velocities  by  observing  the  travel  times  of  body  waves  which  emerged 
directly  from  the  detonation.  Following  this,  a  reflection  survey  will 
be  made  with  the  vibrator  source  at  the  surface  and  geophones  placed  in 
drill  holes  which  extend  down  to  the  altered  layer.  The  vibrator  will 
be  run  from  both  ends  of  this  geophone  spread.  A  survey  of  this  type 
will  be  especially  useful  as  a  check  on  validity  of  the  tunnel  response 
model  since  surface  waves  and  the  strong  P  refraction  arrival  will  be 
attenuated. 

The  geophones  will  then  be  removed  from  the  holes  and 
placed  on  the  surface.  The  spacing  of  a  cluster  of  geophones  at  the 
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surface  will  be  optimized  to  attenuate  Rayleigh  waves.  The  vibrator 
will  again  be  run  from  both  ends  of  the  line.  Paper  and  magnetic  tape 
records  will  be  obtained  from  these  tests.  Evaluation  of  this  reflec¬ 
tion  survey  to  optimize  the  holography  experiment  will  include  travel 
time  corrections,  attenuation,  time  gating  criteria  for  pulse  mode  ex¬ 
periments,  and  stacking  requirements. 

4.5  EQUIPMENT  REQUIREMENTS  FOR  HOLOGRAPHIC  FIELD  EXPERIMENT 

Essentially  all  the  equipment  and  instrumentation  required  to  per¬ 
form  a  seismic  holography  field  experiment  is  presently  available  at 
United  Geophysical  Corporation,  a  subsidiary  of  The  Bendix  Corporation. 
Modifications  of  an  existing  vibrator  source  (described  in  Section  5.1) 
is  required.  The  design  and  fabrication  of  phase  detection  and  integra¬ 
tion  instruments  for  the  recording  of  CW  holographic  data  is  also  re¬ 
quired.  Similar  instrumentation  har  been  designed  and  used  by  the  Bendix 
Research  Laboratories  in  the  past,  however,  so  a  minimum  of  design  effort 
is  required  here. 

A  brief  description  of  some  of  the  key  equipment  available  for 
this  experiment  is  given  below. 

4.5.1  Vibrator  Source  for  Holographic  Surveys 

The  vibrator  system  which  will  be  used  for  all  holographic 
surveys  is  a  modified  version  of  an  existing  unit  operated  by  United 
Geophysical  Corporation.  The  unmodified  and  modified  vibrator  character¬ 
istics  are  given  in  Section  5.1. 

The  modified  vibrator  system  will  provide  an  actuator  out¬ 
put  of  full  rated  force  (10,000  pound  vector)  at  50()  Hz.  The  displace¬ 
ment  amplitude  of  the  reaction  mass  will  be  5  x  10"  inch.  At  200  Hz 
the  expected  amplitude  is  2.5  x  10'  ^  inch  and  at  100  Hz,  1  x  10”^  inch. 

The  vibrator  can  give  a  continuous  wave  (CW)  output,  pulse 
modulated  (PM)  or  sweep  frequency  (chirp)  output  depending  upon  survey 
needs . 


4.5.2  Geophones  for  Holographic  Surveys 

Geophones  of  the  required  capability  (i.e.,  high  voltage 
output,  and  sensitive  to  high  frequencies)  are  readily  available.  The 
Geo  Space  Corporation  Type  GSC-11D  is  an  example  as  is  the  Electro- 
Technical  Labs.  Division  EVS-2  detector.  For  placement  in  drill  holes, 
special  cases  are  available  to  protect  the  geophone  mechanism  from  water 
intrusion.  For  placement  at  the  surface,  a  normal  Spike  base  will  be 
used.  The  United  Geophysical  Corporation  uses  such  geophones  in  its 
daily  operations. 
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•  Characteristics  of  the  GSC-11D 

Natural  undamped  frequency,  10  Hz 

Frequency  response  essentially  flat  from  30  Hz  tc  500  Hz 
Sensitivity,  0.6  V/in/sec 


Characteristics  of  Electro-Tech  EVS-2 
Natural  undamped  frequency,  20  Hz 

Frequency  response  essentially  flat  from  100  Hz  to  300  Hz 
Sensitivity,  0.4  V/in/sec 


4,5,3  Amplifier  and  Recording  System  for  Holographic  Surveys 

Recording  seismic  signals  will  be  carried  out  using  an 
amplifier  system  of  the  United  Geophysical  Corporation.  The  unit  is  a 
Geo  Space  Model  211  High  Gain/Low  Distortion  Digital-Analog  Seismic 
Amplifier  System.  This  system  combines  24  channels  of  seismic  auxiliary 
channels,  control  circuitry  including  gain  programming.  Solid  state  com¬ 
ponents  and  circuit  techniques  are  used.  The  amplifier  fully  meets  the 
needs  of  digital  recording  and  analog  recording.  The  system  can  be 
operated  with  automatic  gain  control  or  in  a  programmed  gain  mode. 

Some  characteristics  of  the  system  are: 

Frequency  Response,  3  to  above  500  Hz 


Equivalent  Input  Noise,  0.1  microvolt 
Dynamic  Range,  120  dB 


M  .  ,  ,,no  J!*6  recordln8  system  is  a  Leech  Control  Recorder  System 
Model  3309.  This  system  includes  a  Model  339  Pulse  Code  Modulation 

AyX«mwnndB1  325°  Mf8netic  Tape  Recorder,  Model  3330  Power  Supply  and 
a  MTR-3200  Recorder-Reproducer  Tape  Deck. 

Some  characteristics  are: 


Sampling  Rate,  1  millisecond 

Dynamic  Range,  78  dB 

Recording  Time,  3  sec  to  10  sec 
in  increments  of  1  sec  or  Manual 


4.6  SUMMARY  OF  EXPERIMENTS 

It  is  recommended  that  an  experimental  field  test  be  oerformed  to 
-:luae  the  capability  of  holographic  techniques  in  detecting  and  imag- 

in8Spr?r8rr^d  °5J/C/t8'  7he  exPerlments  to  be  performed  are  discussed 
in  Section  4.3  and  4.4  and  include  preliminary  refraction  and  reflec- 

„  °nJXP!r  ment8,Jt°  °Ptlml2e  the  hol°8raPh±c  ^ata  acquisition.  Holo- 
g  phic  data  would  be  obtained  using  monofrequency  continuous  wave 
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excitation,  pulsed  monofrequency  excitation,  frequency  modulated  (chirp) 
excitation,  and  quasi-monofrequency  excitation  using  impulsive  (explosive) 
sources.  A  comparison  of  the  results  using  these  different  excitations 
should  give  an  insight  into  the  merits  of  each. 

It  is  estimated  that  a  field  experiment  of  this  magnitude  will  re¬ 
quire  about  one  month  to  perform,  total  time  required  to  prepare  for 

the  experiment  and  to  process  and  analyze  the  data  would  be  aproximately 
one  calendar  year, 

A  possible  daily  schedule  of  activities  during  the  field  experiment 
is  given  below. 


Day 

1. 


2. 


3. 


Operations 

a.  Move  on  site. 

b.  Survey  Hard  Hat  refraction  line. 

c.  Test  vibrator  and  recording  truck  systems. 

a.  Refraction  survey  of  Hard  Hat  site. 

b.  Initial  geophone  spread  based  upon  previously 
established  geologic  model. 

c.  Survey  Piledriver  refraction  line. 

a.  Drill  Hard  Hat  reflection  line.  Spacing  of 
drill  holes  and  depth  of  drilling  is  deter¬ 
mined  from  geologic  model  and  interpretation 
of  Hard  Hat  refraction  survey. 

b.  Refraction  survey  Piledriver  line. 


a.  Complete  drilling  of  Hard  Hat  reflection  line. 

b.  Emplace  geophones  and  load  shot  holes. 

c.  Start  survey  of  two-dimensional  holographic  array. 

d.  Interpretation  of  Piledriver  refraction  survey  data. 

a.  Hard  Hat  reflection  survey  Part  1.  Shot  in  altered 
layer  and  geophones  at  surface.  Shoot  both  ends 

of  the  line. 

b.  Drill  Piledriver  reflection  line. 

c.  Continue  survey  of  holographic  array. 

a.  Continue  reflection  survey  at  Ha^d  Hat  site 
££££_!•  Vibrator  on  surface  and  geophones 
in  drill  holes. 
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b.  Complete  drilling  of  Piledriver  reflection  line. 
Emplace  geophones  and  load  shot  holes. 

c.  Continue  survey  of  holographic  array. 


Continue  reflection  survey  at  Hard  Hat  site 
Part  3.  Vibrator  on  surface  and  surface  geophones. 
The  number  of  geophones  per  amplifier  and  spacing 
of  geophones  will  be  determined  by  the  geologic 
model  and  modifications  to  that  model  obtained 
from  the  reflection  and  refraction  data. 

Continue  survey  of  holographic  array. 

Continue  reflection  survey  at  Hard  Hat. 

Continue  survey  of  array. 

Reflection  survey  at  Piledriver.  Part  1. 

Shot  in  altered  layer  and  geophones  at  surface. 

Continue  survey  to  delineate  nuclear  detonation  area 

Reflection  survey  at  Piledriver.  Part  2. 

Vibrator  at  surface  and  geophones  in  altered  layer. 

Continue  survey  of  holographic  array. 

Reflection  survey  at  Piledriver.  Part  3. 

Vibrator  on  surface  and  geophones  on  surface. 

Begin  holographic  data  acquisition.  For  a  square 
array  containing  90C  sampling  points,  by  sampling 
24  points  at  a  time  (one  line  of  geophones) , 
approximately  9  days  will  be  required  to  sample 
the  entire  aperture.  This  allows  about  4  moves 
of  the  geophone  line  per  day. 

Move  out  with  equipment. 
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SECTIOi,  5 

EVALUATION  OF  SEISMIC  ENERGY  SOURCES  AND  DETECTORS 


Presently  available  commercial  seismic  energy  sources  and  detectors 
were  evaluated  to  determine  if  their  characteristics  were  consister.c  with 
system  requirements.  The  characteristics  of  the  useful  available  or 
modified  components  which  were  used  in  the  design  of  the  holographic 
experiment  were  discussed  in  Section  4. 

5.1  SEISMIC  ENERGY  SOURCES 

Hydraulically  driven  vibratory  sources  have  been  used  in  the  seismic 
exploration  field  for  several  years.  Since  they  are  primarily  used  for 
the  petroleum  industry  and  the  need  Is  to  define  the  earth's  layering 
structure  at  great  depths,  these  sources  are  generally  designed  to  oper¬ 
ate  at  frequencies  from  about  5  to  less  than  100  Hz. 

Two  approaches  can  be  taken  to  obtaining  higher  frequency  sources 
capable  of  generating  a  monofrequency  signal:  a  new  type  of  source 
might  be  designed,  or  presently  available  vibrators  might  be  modified 
making  them  capable  of  higher  frequency  operation.  Since  the  first 
approach  would  be  expected  to  more  costly,  an  effort  was  made  to  see 
if  the  second  approach  was  teasible. 

It  was  determined  that  presently  available  vibrators  could  be 
modified  to  operate  as  high  as  500  Hz.  Shore  Western  Co.,  Monrovia, 
California,  has  had  experience  in  modifying  vibrator  assemblies  to  oper¬ 
ate  at  higher  frequencies.  Discussions  with  Shore  Western  resulted  in 
their  proposing  the  modification  of  a  Wabco  Vibrator  Model  600  B-D  (used 
by  a  subsidiary  of  the  Bendix  Corporation)  capable  of  operating  in  the 
frequency  range  from  5  to  500  Hz. 

The  unmodified  vibrator  is  described  below. 

Vibrator  Model  600  B-D  mounted  on  a  special  off-road  vehicle. 

Electronics:  (Solid  State) 

Designed  by  Bendix  Research  Laboratories 
Phase  Compensator  -  Zero  Crossover  Principle 

Wabco  Vibrator  Assembly: 

Frequency  range  to  200  Hz 

Servo  Valve  -  Hydraulic  Controls  DS4-50G 
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Weight  Assembly  -  Wabco  Low  Profile; 

5  Square  Inches  Ram  Area; 

Vertical  Displacement  +1.25  inch 

Reactionary  Mass  -  4000  pounds 

Base  Plate  -  83.5  inches  x  36.75  inches 
Air  Bag  Isolation 

Base  Plate  Hold  Down  Weight  -  2  1000  pounds 
Hydraulic  Power  Supply 
GMC  Diesel  6-21 

Kline  Pump  Model  1265,  110  GPM 

The  modifications  required  to  provide  5-500  Hz  operating  capabilities 
are  listed  below. 

(1)  Installation  of  a  manifold  and  high  frequency  servovalve 
vMode 1  SV-220)  in  lieu  of  the  existing  valve 

(2)  Increase  Ram  area  to  9  square  inches 

(3)  Install  an  additional  electronic  amplifier  (200  Watt  peak) 
to  drive  the  new  valve 

(4)  Implement  a  second  feedback  loop  based  on  slave  valve 
(second  stage)  position 

(5)  Perfoini  major  modifications  to  the  ground  contact  pad  to 
transmit  the  500  Hz  vibration  into  the  ground 

(6)  Reinforce  vibrator  structure  to  assure  that  high  frequency 
vibration  will  not  cause  failures  when  subjected  to 
extended  field  use 

(7)  Establish  the  exact  performance  of  the  unit  in  terms  of 
ground  vibration  versus  frequency 

These  modifications  to  the  system  will  provide  an  actuator  output 
of  full-rated  force  (10000  pounds  vector)  at  500  Hz.  The  displacement 
of  the  reaction  mass  will  be  5  x  ]0“5  inch.  At  200  Hz,  the  expected 
displacement  amplitude  is  2.5  x  10”3  inch  and  at  100  Hz  1  x  10“ 2  inch. 

5.2  SEISMIC  DETECTORS  -  GEOPHONES 

Commercially  available  geophones  used  in  the  seismic  exploration 
field  wete  found  to  be  adequate  for  the  holographic  system  requirements. 
Since  the  minimum  useful  output  of  geophones  is  related  to  the  noise 
level  of  seismic  amplifiers  (typical  amplifier  noise  is  0.1  pV) ,  geo¬ 
phones  capable  of  detecting  ground  displacements  resulting  in  an 
0.1  pV  output  are  sufficient. 
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One  such  geophone  is  the  Geo-Space  Type  GSC-1IU  detector.  It  is 
a  velocity-sensitive  phone  and  has  a  fiat  response  o*  0.6  tf/in/sec  over 
the  frequency  range  20  to  500  Hz. 

Since  a  minimum  output  of  0.1  uV  is  required,  the  minimum  velocity 
must  not  be  less  than 


v 


io-7  V 


0.6  V/in/sec 


£  1.7  x  10 


-7 


in/sec 


Since  v  =  2  n  f  d ,  where  d  is  ground  displacement,  the  minimum  displacement 
required  to  obtain  0.1  uV  is 

d(100  Hz)  =  2.7  x  10"10  inch 

d(200  Hz)  =  1.35  x  10_1°  inch 

d(500  Hz)  =  5.4  x  10_U  inch 


These  minimum  displacements  indicate  that  both  tunnels  at  the  field  test 
site  (see  Table  4-1)  are  not  detectable  at  500  Hz,  but  both  can  be 
detected  at  lower  frequencies. 
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SECTION  6 

HOLOGRAPHIC  WEAK  SIGNAL  ENHANCEMENT  STUDIES 


The  development  of  a  holographic  weak  signal  enhancement  technique 
was  initiated  at  Bendix  Research  Laboratories  prior  to  the  award  of  this 
contract.  Further  studies  of  this  technique  as  it  pertains  to  seismic 
hologrqnhy  were  made  during  this  program.  This  technique,  which  lends 
itself  wO  computer  simulation,  has  shown  promise  in  defining  images 
of  targets  that  have  signal  strengths  too  weak  to  be  seen  by  conventional 
holographic  imaging  (reconstruction)  techniques,  as  is  often  the  case 
when  a  scattering  target  is  located  above  or  below  a  strong  signal 
reflecting  layer. 

The  theory  used  to  develop  this  technique  and  its  implementation 
were  described  in  a  paper  by  R.  K.  Mueller,  R.  R.  Gupta,  and  P.  N.  Keating, 
"Holographic  Weak-Signal  Enhancement  Technique,"  J.  Appl.  Phys .  ,  Vol.  43, 
February  1972.'  7  Rather  than  duplicating  much  of  this  paper  here,  a  copy 
is  included  in  this  report  as  Appendix  A. 

During  this  program,  two  different  computer  models  representing  a 
two-layer  model  of  the  earth  were  investigated.  The  first  model  (Figure 
6-1)  contained  a  weak  scattering  target  below  the  reflecting  interface 
of  the  two  layers.  This  case  closely  corresponds  to  the  earth  model 
developed  in  Section  2  to  represent  the  selected  field  test  site.  The 
second  model  contained  a  weak  scattering  target  above  the  reflecting 
layer.  In  the  computer  simulations,  all  mode  conversions  (for  example, 
P-waves  to  S-waves)  are  allowed  at  both  the  reflecting  and  scattering 
surfaces  except  when  otherwise  specified.  Lossless  media,  however,  are 
assumed  for  both  models. 

The  scattered  fields  in  the  holographic  plane  are  based  upon  the 
theoretical  analysis  of  Knopoff^  who  has  derived  expressions  for  the 
scattering  coefficients  for  a  plane  P-wave  incident  on  a  perfectly  rigid 
spherical  obstacle.  The  sphere  is  assumed  small  compared  to  the  incident 
wavelength.  Knopoff's  analysis  finds  that  the  scattered  P-wave  spatial 
distribution  is  doubly  circular,  equal  in  the  forward  and  back  directions 
with  nulls  normal  to  the  direction  of  the  propagation.  The  scattered  S- 
wave  distribution  is  also  doubly  circular,  with  a  maximum  amplitude  about 
three  times  greater  than  the  scattered  P-waves,  but  rotated  in  space 
through  90  uegrees  so  the  null  lies  in  the  direction  of  propagation. 

Results  using  the  first  model  are  shown  in  Figures  6-2  through 
6-4.  Figure  6-2  shows  a  conventional  reconstruction  of  the  strong  reflect¬ 
ing  target  (hereafter  called  S\)  and  the  weak  scattering  target  (here¬ 
after  called  S2)  in  the  image  plane  of  S^.  Only  the  image  of  Sj[  can  be 
seen.  Figure  6-3  is  again  a  conventional  reconstruction,  this  time  in 
the  image  plane  of  S2 .  Again  only  the  out-of-focus  image  of  is 


6-1 


Figure  6-1  -  Layout  in  the  Computer  Simulation  Model 


visible.  Figure  6-4  shows  the  results  when  weak  signal  enhancement  tech¬ 
niques  (WSET)  are  applied.  The  image  of  S2  is  clearly  visible,  along 
with  its  conjugate  image,  and  the  image  of  Si  is  suppressed. 

The  results  using  the  second  model  ^Figure  6—5)  are  shown  in 
Figures  6-6  through  6-15.  Figure  6-7  shows  a  conventional  reconstruction 
for  the  case  of  a  single  scattering  target.  The  reconstruction  is  in 
the  plane  of  the  scatterer.  Only  the  out-of-focus  image  due  to  the 
reflecting  layer  can  be  seen;  S2  cannot  be  resolved.  Figure  6-7  shows 
the  result  for  the  same  case  when  WSET  is  used  and  only  scattered  P-waves 
are  allowed.  The  image  S2  is  now  well  resolved.  Figure  6-8  is  the  result 
for  the  same  case  as  Figure  6-7  except  that  random  scattering  targets 
(noise)  are  now  included  in  the  model.  Again,  good  image  resolution 
is  obtained.  Figure  6-9  again  gives  the  results  for  a  single  scatterer 
(S2) ,  but  here  both  scattered  P-  and  S-waves  are  allowed.  Jn  this  case 
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tMHANCED IMAGE 


6-5 


S2  is  Much  Larger  Than  Sj  After  Enhancement. 
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ONLY  SCATTERED  P-WAVES  ARE  ALLOWED. 
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Figure  6-7  -  Weak  Signal  Enhancement  Reconstruction  of  a  Single  Scatterer 
in  Image  Plane  of  S2*  Only  Scattered  P-Waves  are  Allowed. 


RANDOM  SCATTERING  TARGETS  ARE  INCLUDED  IN  THE  MODEL. 
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Figure  6-8  -  Weak  Signal  Enhancement  Reconstruction  of  a  Single  Scatterer 
Image  Plane  of  S2-  Random  Scattering  Targets  are  Included  in  the  Model. 


BOTH  SCATTERED  P-  AND  S-WAVES  ARE  ALLOWED. 
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interference  between  S-  and  P-waves  greatly  reduce  the  resolution  of  the 
reconstructed  image . 


When  an  extended  scatterer  (S2)  is  used  instead  of  a  single  scatterer, 
the  results  are  as  follows.  The  conventional  reconstruction  is  essen 
tially  the  same  as  shown  in  Figure  6-7  for  a  single  scatterer,  i.e.,  2 
cannot  be  resolved.  Figure  6-10  shows  the  results  when  only  scattere 
P-waves  are  allowed  and  with  no  random  scatterers.  Figure  6-11  is  tor  tne 
same  case,  except  random  scatterers  are  included  in  the  model.  n 
cases,  S2  is  again  well  resolved.  Figure  6-12  gives  the  results  when 
both  scattered  P-  and  S-waves  are  allowed.  The  S-wave  interference  i 
again  apparent  in  that  it  greatly  reduces  the  S2  image  resolution. 

The  results  given  above  indicate  that,  at  least  for  the  model  shown, 
interference  due  to  scattered  S-waves  influence  the  resolution  of  the 
desired  image  to  a  greater  extent  than  does  the  presence  of  ran  om >  nout 
sources.  The  computer  model  used,  in  which  a  lossless  medium  is  assumed, 
is  actually  a  worse  case  condition  with  respect  to  S-wave  interference. 

In  a  lossy  medium,  S-waves  typically  attenuate  more  than  P-waves,  wnich 
should  result  in  less  S-wave  interference  at  the  hologram  plane.  Examples 
of  the  increase  of  resolution  obtained  when  shear  waves  are  allowed  to 
attenuate  more  than  congressional  waves  are  shown  in  Figures  6-13  th.ough 
6-15.  In  Figure  6-13  the  shear  waves  are  attenuated  by  a  factor  ct  , 
in  Figure  6-14  by  a  factor  of  4,  and  in  Figure  6-15  by  a  factor  oi  8. 

As  would  be  expected,  the  image  resolution  improves  with  shear  wave 

suppression . 


The  results  given  in  this  section  indicate  that  the  weak  signal 
enhancement  technique  can  be  a  powerful  tool  in  the  reconstruction  of 
seismic  holograms. 
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ONLY  SCATTERED  F-WAVES  ARE  ALLOWED. 


I  tlOl  IY4 


in 


Figure  6-10  -  Weak  Signal  Enhancement  Reconstruction  of  Extended  Scatterers 
in  the  Image  Plane  of  S2.  Only  Scattered  P- Waves  are  Allowed. 
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BOTH  SCATTERED  P-WAVES  AND  S 
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SCATTERED  S-WAVE  ATTENUATED  BY  FACTOR  OF  2 
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in  Image  Plane  of  S2-  Scattered  S -Wave  Attenuated  by  a  Factor  of 


SCATTERED  S-WAVES  ATTENUATED  BY  FACTOR  OF  4. 
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Figure  6-14  -  Weak  Signal  Enhancement  Reconstruction  of  an  Extended  Scatterer 
in  Image  Plane  of  S2.  Scattered  S-Waves  Attenuated  by  a  Factor  of  4. 


SCATTERED  S-WAVES  ATTENUATED  BY  FACTOR  OF  8. 
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Figure  6-15  -  Weak  Signal  Enhancement  Reconstruction  of  an  Extended  Scatterer 
in  Image  Plane  of  S?.  Scattered  S-Waves  Attenuated  by  a  Factor  of  8. 
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APPENDIX  A 

HOLOGRAPHIC  WEAK  SIGNAL 
ENHANCEMENT  TECHNIQUE 


Holographic  Weak-Signal  Enhancement  Technique 


R.  K,  Mueller,  R.  R.  Gupta,  and  P,  N.  Keating 
Bendix  Research  Laboratories,  Southfield ,  Michigan  48076 
{Received  t6  June  1971) 

A  holographic  weak-aignal  enhancement  technique  {WSET)  ia  described  hich  enhaneea  the 
image  of  the  source  of  a  weak  algnnl  of  intercat  in  the  preaenee  of  an  unwanted  atrong  signal 
with  very  ltttle  a  priori  knowledge  of  the  latter.  The  technique  ia  most  auitable  for  imaging 
applications  in  fields  such  as  aeiamic  or  oceanographic  holography,  where  the  desired  weak- 
scattered  signal  is  often  overwhelmed  hy  a  atrnng  reflected  signal  from  interfaces  between 
different  media  The  feasibility  of  the  WSET  has  been  shown  by  simulating  a  two-sol  id- lay¬ 
er  acoustic  problem  on  a  digital  computer.  Well-defined  enhanced  Imagea  have  heen  ohtained 
for  target  scattering  which  ia  ao  much  weaker  thnn  a  strong  reflected  signal  that  the  target 
imagea  were  lost  in  noise  using  conventional  holographic  reconstruction. 


I  iN  f  ROUUITtON 

imaging  systems,  as  all  other  information-carrying 
systems,  are  plagued  with  the  problem  of  identifying 
weak  signals  in  the  presence  of  noise.  For  example, 
"noise”  front  the  wings  cf  a  nearby  strong  signal  often 
dominates  and  obliterates  a  weak  signal,  especially  if 
the  former  is  out  of  focus  h.  the  desired  Image  plane. 
This  type  of  weak-signal  Identification  problem  fre¬ 
quently  arises  when  conventional  holographic  Imaging 
techniques  are  used  to  image  utderground  or  ocean¬ 
ographic  weak-scaiierlng  targets.  In  these  applications 
the  weak  desired  signal  scattered  from  tie  iargets  can 
be  lost  In  the  background  of  the  strong  signal  due  io  the 
reflection  of  the  propogating  waves  from,  for  example, 
the  stratified  layers  of  the  earth,  or  from  ihe  ocean 
bottom.  In  addition,  there  may  be  noise  due  to  other 
sources,  such  as  propagating  surface  waves.  In  such 
applications  It  Is  very  desirable  to  have  a  processing 
technique  which  could  enhance  the  (desired)  weak  Image 
in  the  presence  of  a  strong  (unwanted)  signal  of  simple 
form.  For  convenience,  the  unwanted  strong  signal  is 
referred  to  as  a  reflected  signal  unless  otherwise 
stated,  throughout  the  article. 

In  this  paper  we  describe  a  new  weak-signal  enhance¬ 
ment  technique  (WSET)  which  enhances  the  Image  of 


weak  iargets  when  accompanied  by  a  strong  reflection, 
with  virtually  no  a  priori  knowledge  of  the  scene.  The 
WSET  would  be  applicable  to  problems  where  the  fol¬ 
lowing  two  condlilons  are  met:  (i)  The  unwanted  signal 
is  cohereni  with  the  target  signal  and  Its  intensity  at  the 
holographic  plane  Is  relatively  slowly  varying  (ii)  Some 
spatial  separation  of  the  weak  target  and  the  strong-sig¬ 
nal  virtual  source  is  obtainable.  The  first  condition  is 
likely  to  be  met  In  a  wide  class  of  cases;  for  example, 
for  ihe  strong  reflection  from  an  interface.  In  the  case 
of  such  a  strong  signal,  the  second  requirement  can  be 
mei  by  placlngihe  receive  aperture  so  that  the  virtual 
source  associated  with  this  reflection  is  suitably  sepa¬ 
rated  from  ihe  target  in  the  plane  parallel  to  the  holo¬ 
graphic  plane 

The  enhancement  method  Is  based  on  (i)  ihe  use  of  the 
signal  field  for  reconstruction  and  (11)  high-pass  filtering 
of  some  of  the  information.  This  reconstruction  method 
Is  most  conveniently  utilized  by  means  of  digital-com¬ 
puter  reconstruction  techniques. 

Section  II  consists  of  a  description  of  the  proposed  en¬ 
hancement  technique  together  with  an  evaluation  of  Us 
usefulness.  In  Sec.  m  we  present  an  account  of  a  com¬ 
puter- simulated  application  of  the  enhancement  tech¬ 
nique  to  a  geophysical  holographic  problem.  Section  IV 
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consists  of  a  discussion  of  the  simulation  results  and  an 
evaluation  of  the  method  for  actual  field  data. 


II  *LAK  SIGNAL  tMlANC t.WtNT  Tt(  HNIQUfc 

The  enhancement  reconstruction  technique  is  an  exten- 
s  ion  u *  me  lu.„ tntionai  holographic  reconstruction 
method.1  It  consists,  in  essence,  of  tuj  successive 
reconstructions,  the  second  using  the  strong  return  as 
the  reconstruction  signal. 

In  a  conventional  hologram  the  image  information  ts 
stored  as  the  product  R*St  where  S  S,  *  St  is  the  total 
signal  field  at  the  hologram  plane.  Sl  is  the  strong  un¬ 
wanted  signal,  S2the  weak  signal  of  interest,  and  R  is 
the  reference  field.  Conventional  reconstruction  in¬ 
volves  an  addition.*!  multiplication  by  R*  ,  a  reconstruc¬ 
tion  field  at  the  mlograni  plane,  followed  by  a  transfor¬ 
mation  {e.  g.  .  a  Courier  transform  in  the  Fraunhofer 
case,  a  Fresnel  transform  in  the  Fresnel  approxima¬ 
tion)  to  the  image  plane.  With  a  suitable  Rr  ,  this 
transformation  yields,  at  thj  image  plane,  the  field 
S  S  ,  .  s2  which  is  the  image  corresponding  to  the 
holographic  field. 

The  weak- signal  enhancement  technique  involves  the 
following  steps:  (a)  recording  of  S,  as  in  conventional 
holography;  (b)  the  formation  of  IS  I  *  -  iS,  +  S2l2from 
S,  or  directly;  (c)  high-pass  filtering  of  I  SI1  (the 
result  is  denoted  as  1 5 .  (d)  multiplication  of  1 5 1  * 

by  S  and  reconstruction.  U  |S,I2  is  slowly  varying,  and 
S,  and  St  are  spatially  separated,  we  can  carry  out 
high-pass  filtering  (step  (c)|  to  obtain  the  reduction 

•5,|  4  •  15,1 2  •  2  Re(S,$2*  )  -  l$2l2*  2Re(S,V>  •  (D 
Hence,  step  (d)  yields 

/  (5 1  «S2)  (|S:|2  .Si  S2*  *  S,  *S2) 

(IS,!2.  I52|2)s2.2  |S2|  2S,  .  Sf  s\.sfs$.  (2) 

The  last  two  terms  give  signals  in  the  image  plane  which 
are  spatiaiiy  separated  from  the  images  of  interest  and 
will  not  be  considered  further.  The  third  term  yields 
essentially  the  image  S,  but  it  is  reduced  by  the  small  fac  - 
tor  2;S2i2.  The  first  term  is  the  desired  image  term 
multiplied  by  the  large  factor  IS,  2  (since  \h :  ■ 2 1 S2 1 2). 
Thus,  the  weak  image  has  been  enhanced  reiauve  to  the 
st  rong  one  by  the  large  factor  J  :S,  i2 '  IS2| 2. 

The  condition  that  I  S,  | 2  is  a  siowly  varying  intensity  in 
the  holographic  plane  is  necessary  so  that  (i )  rejection 
of  |S,I  1  by  high-pass  filtering  is  effective;  (ii)  the  final 
enhanced  fmage  S2  Is  not  degraded.  The  condition  that 
,  i2  are  spatially  separated  is  necessary  so  that 
S,S2  is  a  term  with  high  spatial  frequencies  and 
can  be  readily  filtered  from  I  Sl 2. 

It  is  instructive  to  consider  the  enhancement  technique 
in  the  case  of  point  objects  and  Fraunhofer  holograms. 

In  this  case,  S,t  S2,  and  R  are  all  sinusoidal  functions, 
with  spatial  frequency  determined  by  the  angular  posi¬ 
tions  q  ••  ((sin  P)/x  |  (cos  d>,  sin  <M  of  the  point  sources, 
where  9  and  d>  are  the  polar  and  azimuthal  angles, 
respectively.  For  example,  S,  -  exp(iq,  x).  S2 
exp(icf2.  x  )(o^  1),  where  x  is  the  two-dimensional 


position  vector  In  the  hologram  plane.  Thus 
1 5| 2  =  f  1  *  o2  ♦  2o  cos(4 1  -  52)  •  x), 

ISI2  2  a  cos(tf  -  q2)  •  x  , 
if  the  dc  is  rejected  andq,*q2.  Hence 
/  2a  f  exptiq,  *x  )  *  o  expf/q  ,•  x)|  cos(q,  -  52). x 

of  exp(iq2  •  x  ) »  o  exp(i  q, .  x  )  »  exp|  i(2<fi  -  q^  •  x  J 
♦  o  exp|  i(2q2  -q,)-  x||  .  (4) 

We  thus  obtain  four  point  images.  The  amplitude  of 
J>2  is  now  larger  than  that  of  s,  by  a  factor  of  1  'a ,  in¬ 
stead  of  being  smaller  by  a  factor  of  o.  The  other  two 
images  occur  at  positions  different  from  s2,  pro¬ 
vided  qx  *  <7:  (i.  e. .  the  images  are  not  coincident),  as 
already  assumed.  (We  note  that  in  this  case  i  S2| 2  is 
also  filtered  out  so  that  the  additional  factor  of  2 
vanishes.  ) 

Steps  (a)- (d)  necessary  for  the  enhancement  technique 
can  be  carried  out  with  optical  components,  although  not 
readily.  They  can.  however,  be  very  easily  cai  ried 
out  on  a  digital  computer,  where  effective  filtering  and 
large  dynamic  range  are  readily  available. 

In  a  problem  where  the  equivalent  source  of  the  signal 
S,  is  a  point,  step  (a)  can  be  eliminated  by  using  a 
synthetic  point-source  signal  (instead  of  S  as  proposed) 
in  step  (d)  lo  illuminate  |S|2.  However,  in  cases  where 
S,  is  not  from  a  point  source,  as  in  the  present  study, 
it  is  preferable  to  retain  step  (a)  and  use  this  to  recon¬ 
struct  the  image, 


111  COM  PUT  tR  SIMULATION  01  V\L  AK  SIGNAL  t  Ml  AMT  MINT 
TLCHNlQUL 

A  Modi* t  amt  Computational  Method 

Since  the  experimental  evaluation  of  the  feasibility  of  the 
proposed  holographic  weak-signal  enhancement  technique 
in  seismic  or  oceanographic  holography  would  be  ex¬ 
pensive  and  time  consuming,  a  simple  (idealized)  case 
of  seismic  holography  was  first  simulated  on  a  digital 
computer.  In  this  simulation,  it  was  desired  to  use  a 
digital  computer  to  reconstruct  the  image  of  a  small 
acoustic  scatterer  buried  below  an  interface  between  two 
acoustic  media  (as  shown  in  Fig.  1);  the  source  and 
linear  detector  array  are  above  the  interface.  The 
scattering  strength  of  the  target  was  chosen  such  that 
the  conventionally  reconstructed  holographic  image 
«S2  is  not  identifiable  in  the  presence  of  the  noise  ac¬ 
companying  the  strong  image  In  the  computer  sim¬ 
ulation,  we  have  used  the  following  model: 

(i)  The  transmitting  (acoustic)  source  radiates  coherent 
longitudinal  waves  (i.  e. ,  P  waves),  (ii)  The  interface 
between  the  two  media  is  planar,  (iii)  Both  media  are 
homogeneous,  isotropic,  and  lossless.  ( i v )  Both  media 
are  solid,  and  hence  mode  transformation  from  P  waves 
to  S  waves  (i.  e.  ,  shear  waves)  and  vice  versa  takes 
place  whenever  propagating  waves  are  either  reflected, 
refracted,  or  scattered,  (v)  Either  no  surface  waves 
reach  the  detector  array  or  they  can  be  discriminated 
against  by,  for  example,  gating,  (vi)  The  holographic 
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fields  are  detected  by  a  finite  linear  array  of  detectors 
located  at  many  wavelengths  away  from  the  interface, 

(vul  The  sc  at!  erer  and  virtual  source  are  located  many 
wavelengths  away  fr>»m  the  interface,  vet  they  are  in  the 
near-field  region  of  the  receive  aperture,  as  would  be 
the  usual  case  in  seismic  holography.  Although  ray- 
optics  approximations  t * »  ihe  solution  of  propagating 
spherical  waves  are  valid  ir.  the  far- field  (and  not  the 
near -field)  region  of  the  aperture, 2  m  our  simulation 
study  we  shall  for  simplicity  use  the  ray -optics  ap¬ 
proximations.  (viii)  There  is  sufficient  angular  sepa¬ 
ration  between  the  sc  at  t  erer  and  the  unwanted  strong- 
signal  s  »urce  so  that  bev -frequency  co moments  can 
be  readily  filtered  "ui  from  the  !.ST  signal. 

Although  the  assumed  two- layer  model  for  the  earth 
is  nui  entirely  realistic,  the  conclusion  arrived  at  here 
will,  howexer.  Le  .ippl  tcahle  tn  a  seisi  uc  proh, cm  where 
the  reflect  cd  signal  S,  ts  caused  predominantly  by  a  single 
interface  We  note  however  that  the  inclusion  of  both  S 
and  P  waxes  results  m  quite  a  stringent  test  of  ilu* 
nielli  *d  since  the  strong  xirtual  source  is  a  compound 
ob|c\  i  and  the  mterlerence  between  the  S  and  /*  reflected 
waxes  must  also  be  tillered  out. 

The  computer  simulation  of  such  a  problem  is  therefore 
expected  n  provide  a  reasonable  evaluate  n  of  the  pro¬ 
posed  holographic  weak-signal  enhancement  technique 
when  applied  to  a  practical  problem. 

The  Rayleigh  - Sommer  fold  diffraction  integral  formula3 
(  ».  (5)]  is  used  to  reconstruct  the  original  scene  from 

tu.  t  Ids  detected  i»v  the  linear  array.  For  the  general 
confihurat i  'll  shown  in  Fig.  2,  the  diffraction  integral 
formula  calculates  the  held  amplitude  '(i,i  . /)  at  the 
jsiint  /Mi ,  v  ,  i  >  from  the  known  field  distribution  /'  ,  n  l 

over  the  finite  planar  aperture 

The  integral  formula  is3 

f\ i.x.^1  (1  /'if  A’trlcusfn- r> Irftrfn  .  (5) 

where 
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Kir)  c  »  (three-dimensional  kernel),  (b) 

r  (i-02-  (y-i| 

>  is  the  propagation  wavelength  of  the  detected  fields, 
and  rnsfii-  ? )  is  the  cosine  of  the  angle  iietw  eon  the  out  - 
ward  normal  n  and  the  vector  r  joining  Q  to  P  This 
formula  assumes  a  scalar  description  of  the  acoustic 
fields.  In  our  simulated  problem  this  would  yield 
reasonably  accurate  results  since  the  signal  at  the  re- 
ceive  aperture  is  mostly  contributed  by  the  longitudinal 
waves. 

In  the  Fresnel  (or  Fraunhofer)  /one,  the  image  may  be 
reconstructed  using  the  simpler  Fresnel  (or  Fraunhofer) 
diffraction  formula,  instead  of  the  Rayleigh  detraction 
formula  which  should  be  used  in  the  near-field  zone.  In 
our  simulation  w*ork  we  are  concerned  with  targets  and 
virtual  sources  which  are  in  the  near  field,  as  is  the 
usual  case  with  the  long  wavelengths  utilized  in  seismic 
holography. 

In  order  to  uniquely  determine  fields  in  an  unbounded 
space  by  using  the  diffraction  integral  formula,  the  field 
distribution  must  be  specified  over  the  two-dimensional 
aperture.  In  our  simulated  problem  we  have,  for 
simplicity,  computed  fields  only  along  a  linear  array. 
Therefore,  for  our  near-field  simulated  problem,  the 
three-dimensional  kernel  A  (W  used  in  diffraction  integral 
Fq  (5)  should  be  modified  to  reconstruct  two-dimen¬ 
sional  Images  from  one-dimensional  holograms  For 
simplicity,  however,  we  have  used  the  unmodified  three- 
dimensi  mal  kernal  | Eq.  (6)].  This  is  expected  to  only 
slightly  impair  the  *' conventional"  aperture-limited 
resolution  of  the  reconstructed  images. 

Fast  Fourler-transforni  techniques,4-6  based  on  the 
Cooley-Tukey  algorithm  (known  as  FKT).  havebeenused 
to  digitally  reconstruct  the  holographic  images  using  the 
diffraction  integral  formula  (5).  The  FFT  is  simply  an 
efficient  method  for  computing  the  discrete  Fourier 
transform  (DFT).  The  FFT  can  be  used  in  place  of  con¬ 
tinuous  Fourier  transform  only  to  the  extent  DFT  can 
be  used,  but  with  a  substantial  reduction  in  computer 
time.  The  FFT  can  be  applied  to  operations  such  as 
computing  a  spectrogram,  convolution,  correlation,  and 
filtering.  In  our  problem  the  FFT  algorithm  is  used  (I) 
to  evaluate  the  integral  equation  (5)  via  the  convolution 
theorem  methods  and  (ii)  to  filter  out  the  low-frequency 
terms  from  a  signal. 
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At  dixcuftM*  ***H#r.  il*  l*t c»m- 

pHintf*  Mte  i » b*  ftltrrwl  »*  Irww  ihr  tt***!  iSl 
jitiwe  m  «ur  fimulitrd  *r  Have  kirntl" 

nfe*  <4  tti#  relative  «4  the  refleeitmt  interface 

ami  veatterer  ih*  tfuiul  tui.4f  fr«*uewy  t»  easily 
determined  THr*r  E<»«.  spinal -frequency  components 
ire  then  Itliered  mui  ln»m  thr  I  SI  *  sttfvul  by  a  dt*ttal 
hi*h-iu*»  filler.  5tme  in  any  practical  imaging  |»rob- 
Irm.  thr  tar*H  l.tcaiMna  arr  nut  known  *  ptitifi  .  ihr 
jppt  »|iriatr  cui«»f|  Ircquency  nuisi  lir  ch-iarn  by  aornr 
altrrnaiivr  mrth««l.  A  Itltrrtn*  approach  that  nay  In* 
fMll.twrd  in  procrsai n«  rral  Itrld  daia  la  presented 
la'rr 

H  Mtnlet  Pjomclff* 

Conventional  holographic  techniques  as  wrll  aa  the 
Ivdographlc  weak-signal  enhancement  technique  arc 
employed  lo  reconstruct  Ihr  original  aernr  ahown  In  the 
configuration  of  Fig.  1  Thr  detectors  in  the  linear 
array  are  assumed  to  be  sensitive  only  to  the  vertical 
component  of  the  force.  The  mathematical  expressions 
lor  the  pressure  waves  resulting  from  either  reflection, 
refraction,  or  scattering  of  waves  are  available  in  the 
lltcrature.2*7'- 

The  following  reasonable  values  of  model  parameters 
have  been  chosen  for  the  evaluation  of  the  proposed 
WSET.  Poisson’s  ratio  for  both  media  -  0.  25:  raiio 
of  densities  of  medium  1  to  medium  II  Is  ratio  of 
velocities  of  P  waves  in  medium  1  to  medium  II  =  1.5} 
ratio  of  the  radius  of  scattering  sphere  to  the  wave¬ 
length*  =  0.  5;  reference  wave  R  -A  exp[(f/v2 )  (v+  a)]; 
length  of  linear  deiccior  array  =64X;  number  of 
detectors  128.  The  size  and  the  locaiion  of  the  deiec- 
tor  array  were  selected  such  that  no  Head-type  surface 
waves2  would  theoretically  arrive  at  ihe  detector  array. 
The  detectors  in  the  linear  are  spaced  ai  an  interval  of 
i*t  which  is  somewhat  less  than  one-half  of  the  small¬ 
est  spatial  wavelengih  X,  in  the  signal  5 ,  so  that  the 
sampling  rate  is  close  to  opiimum.  Wiih  ihe  chosen 
parameter  values,  the  scattered  signal  intensity  \St\ 
was  found  to  be  1.45*  10’4  smaller  than  the  reflected 
signal  Intensity  I  S,| 2. 

C  Dwifn  Performance 

The  image  resoluiion  obtained  is  predominately  aperture 
limited.  [In  addiiion  (as  discussed  earlier),  ihe  use 
of  the  unmodified  three-dimensional  kernel  in  our  cal¬ 
culations  is  expected  to  furiher  impair  the  resoluiion.  ] 

In  our  example  ihe  Images  are  reconstructed  in  ihe 
near-fleld  region  of  ihe  aperture,  where  the  resoluiion 
calculations  are  difficult.  Reasonable  estimaies  for  the 
image  resolution  can,  however,  be  obiained  from  far- 
field  calculations"  if  ihe  Images  are  not  reconsiructed 
quite  close  to  the  aperture,  as  is  ihe  case  in  our  prob¬ 
lem. 

Based  on  far-field  calculations,  ihe  resolvable  linear 
dlsiance  AL  and  the  lateral  distance  (depth  of  focus) 
between  ihe  iwo  point- source  objects  for  an  aperture 
of  width  2a  at  a  distance  z  (normal  to  the  aperiure)  are 

AL  *  X  U/n), 


A/  »  I  Ufmt ,  m 

In  tke process  of  tltfiuJ  enhancement  we  «Maw  ihe  im*** 
<4  thr  Ur get  from  thr  mtsfulsted  aiitial  !.%»•*§  ** 

thr  tig  ail  If  lS|  I1  is  m4  osulint  uvrr  ihr  ht4t«rsphir 
plane,  an  additional  broadening  ts  ihu*  introduced  into 
Ihe  enhancitl  image  However,  since  |S|I*  ts  a  low 
frisjuency  fund  ton.  the  addti  tonal  bias  of  resolution  Is 
not  expected  to  be  !«*>  severe.  Ustng  Fq  (71.  the  ex¬ 
pected  hall-widths  of  it,  and  ihe  enhanced  image  of  *, 
are  rtmghly  J*  .  |* .  resi»ec lively;  Ihe  dcplh  of  focus  of 
the  enhanced  image  is  almul  . 

In  our  cast*,  a  lirst-nrder  estimate  *>f  the  Image  degra¬ 
dation  Introduced  liy  the  usage  of  an  unmodified  kernel 
has  Indicated  a  rather  insignificant  contribution  lo  the 
existing  aperture* limited  resolution  of  the  image 


IV  KtSlHT 
A  Simulated  hiwfilr 

Figure  3  shows  the  imuge  Intensity  distribution  obtained 
by  a  conventional  holographic  reconstruction  of  the  total 
scene  shown  In  Fig.  1.  This  reconstruction  was  done  in 
the  cxpccicd  image  plane  of  the  virtual  source  8t.  The 
Image  of  82  Is  not  Identifiable  from  the  background  noise. 
The  image  location  and  the  resolution  of  the  image  of 
8,  are  as  expected}  the  ratio  of  the  intensity  of  ihe  weak 
signal  to  the  strong  signal  is  about  1.45xlO‘V 

Figure  4  shows  the  image  intensity  obtained  by  another 
conventional  reconstruction  of  the  total  scene  — but  now 
in  ihe  expected  Image  plane  of  s2.  The  expected  Image 
plane  of  s2  is  calculated  to  be  at  a  distance  of  125X  from 
the  hologram  plane.  In  this  plane  the  target  is  located 
ai  a  distance  of  25X  away  from  the  cenier  of  ihe  hologram. 
The  image  from  \SV\*  is  oui  of  focus  and  ihe  image  from 
|  s2|  2  Is  still  unidentifiable  from  ihe  background.  In 
Fig.  5  we  preseni  an  enlargement  of  part  of  Fig.  4  in 
the  vicinity  of  the  expected  locaiion  of  v2.  There  is  no 
trace  of  the  desired  s2  image. 

Figure  6  shows  ihe  weak-signal  enhancement  reconsiruc- 
tion  in  the  expected  image  plane  of  *2  as  would  normally 
be  of  interest.  Again  iniensitles  are  plotted.  The  fig¬ 
ure  shows  (i)  the  image  of  the  weak  target  is  ai  the 
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FIG.  3.  Conventional  holographic  reconstruction  ISj*S2l  in 
the  image  plane  of  S,.  The  desired  image  I  S2 1  2  is  not  obser¬ 
vable  in  the  presence  of  |  St  1 2. 
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expected  locatton,  25*.  and  considerably  enhanced  over 
the  image  of  slf  (ii)  the  half-width  of  the  enhanced  intake 
is  about  ix  (almost  the  same  as  theoretically  predicted) 
and  (ill)  the  weak-target  intensity  ts  enhanced  by  a  fac-  * 
tor  of  S,  I*  I St |4,  The  theoretical  Improvement  obtain¬ 
able  in  the  image  intensities  is  expected  to  be  on  the 
order  of  5>  107.  in  actual  fact,  however,  the  effective¬ 
ness  of  the  filtering  of  \S{\Z  is  likely  to  be  the  limiting 
process. 

Comparison  of  Figs.  3  and  6  shows  the  considerable 
effectiveness  of  the  weak-signal  enhancement  technique. 

It  alters  the  situation  front  one  in  which  the  desired 
I  s al  is  unobservable  in  the  presence  of  J  st  | 2  (Fig.  3) 
to  one  in  which  the  unwanted  intensity  from  s,  is  un¬ 
observable  In  the  presence  of  the  enhanced  image  of  the 
s2  target  (Fig.  6). 

0  Discussion 

The  results  of  the  simulation  study  show  that  targets 


TARGET  COORDINATES  -  IN  WAVELENGTHS 


FIG.  5.  Enlargement  of  part  of  Ftg.  4  in  the  vtctntty  of  the 
expected  iocatton  of  tmage  S2.  1S2|2  ts  attll  undetectable. 


which  trailer  c*>  weakly  to  be  imaged  in  the  presence 
•>f  a  at  rung  reflected  signal  by  conventional  hob  graphic 
methods  can  often  be  readily  imaged  by  the  proposed 
holographic  weak -signal  enhancement  techntqur  The 
enhancement  technique  has  been  found  to  lie  able  to 
image  weak  retuma  that  are  two  to  three  orders  of 
magnitude  below  the  amplttude  of  those  which  are  just 
discernible  by  conventional  holographs  techniques 
A  comparison  of  Ftgs  3  and  6  shows  the  considerable 
enhancement  which  can  be  obtained  (a  factor  of  over 
10'  In  Intensities). 

The  effectiveness  of  the  WSET  depends  upon  the  effici¬ 
ency  with  which  the  high-piss  filter  can  suppress  the 
low-frequency  part  of  |.V| *  The  spatial  frequency  of 
the  filtered  signal  Increase  with  increasing  angular 
separation  between  the  weak  Image  st  and  strong  un¬ 
wanted  Image  Thus  the  effectiveness  of  the  filter, 
and  hence  the  WSET,  improves  with  increasing  angular 
separation  between  the  sources  of  the  desired  weak 
signals,  and  the  undesired  strong  signals.  For  ex¬ 
ample,  when  Sj  is  due  to  a  strong  reflection  from  a 
horizontal  plane,  the  WSET  is  most  effective  when  the 
receivers  are  not  directly  above  the  target9. 

C  .  ProcTMing  of  Field  Data 

All  processing  steps  involved  In  Implementing  the 
holographic  WSET  to  field  data  are  straightforwaid, 
except  for  the  filtering  of  the  low-frequency  components. 
In  the  simulated  case  the  filtering  step  was  readily  car¬ 
ried  out  because  of  a  priori  knowledge  of  the  location  of 
the  reflecting  Interface  In  any  practical  problem,  the 
filtering  step  would,  however,  be  somewhat  more  diffi¬ 
cult,  since  the  locations  of  reflecting  Interfaces  will,  in 
general,  be  unknown  beforehand.  In  a  practical  Imag¬ 
ing  problem  the  following  filtering  procedure  may  be 
followed  lo  Implement  the  WSET. 

(i)  Conventionally  reconstructed  holographic  Images 
are  generated  In  a  suitable  number  of  adjacent 
Image  planes  parallel  to  the  holographic  plane.  These 
Image  planes  should  be  spaced  at  roughly  the  depth  of 
focus  of  the  receive  aperture.  From  these  reconstruc¬ 
ted  images,  the  locations  of  the  focused  strong  images 
rS»  are  obtained,  for  comparison  with  the  enhanced  im¬ 
age  fields  obtained  in  step  (III). 


TARGET  COORDINATES  IN  WAVE  LENGTHS 

FIG.  6.  Weak-signal  enhancement  reconstruction  in  the  tmagi 
plane  of  S2.  After  enhancement,  the  desired  image  I  S,  1 3  has 
become  much  larger  than  I  St  1 3.  (The  same  model  parameters 
have  been  uaed  here  as  for  Ftgs,  3-5.) 
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(ttl  Several  ftllersd  atgnala  tSI*.  each  with  a  dtffrre* 
filter  cutoff  frequency,  are  it  merit  ed 

(ttt)  Emm  each  filtered  atonal  ISI1,  enhanced  tmanea 
of  weak  target*  are  recon*lrucied  in  the  several  image 
planes  by  ustng  the  signal  «,  •  *i  for  reconstruct  ton 
From  the  images  so  reconstructed,  the  optimum  target 
reconstruction*  are  choaen. 

(OM  UNIONS 

We  have  presented  a  technique  which  can  enhance  the 
image  of  a  weak  target  in  the  presence  of  a  strong  (e.g. , 
reflected)  signal,  ily  computer  simulation  we  have 
shown  that  the  proposed  holographic  enhancement  tech¬ 
nique  can  indeed  provide  considerable  enhancement  of 
the  image  of  a  weak  signal  accompanied  by  an  undesired 
strong  signal  coherent  with  the  target  return  signal. 

The  WSET  is  feasible  whenever  (i)  the  strong-signai 
intensity  is  relatively  slowly  varying  in  the  hologram 
plane  and  (ii)  a  sufficient  angular  separation  exists  be¬ 
tween  ihe  target  and  the  virtual  source  of  the  unwanted 
strong  signal,  so  that  effective  filtering  of  the  iow- 
frequency  terms  can  be  achieved.  To  a  great  extent 
the  effectiveness  of  this  filter  determines  the  weakest 
signal  that  can  be  imaged  in  this  manner  by  the  proposed 
WSET.  Based  on  our  computer  simulation  study  we 
calculate  that  in  a  practical  problem  the  proposed  holo¬ 
graphic  WSET  should  readily  image  weak-scattering 


targets  that  would  otherwise  be  lost  in  ihe  background 
of  ihe  image  of  ihe  accompanying  undestred  strong 
coherent  stgnal 

MKNOfclllN.MIM* 

The  author*  are  indebted  to  R  F  Steinberg  and  T 
Sawatari  for  valuable  discuasions  and  J,  Lindsay  tor 
numerical  computations, 


'.I.  W  (HMNtdian  tut riN/*r tton  fo  Fourier  »V>f|rs  IMH I row-ttlll, 
New  Yoik  U«*“>  fhap  s. 

?t.  M  |*n  hhovsklkh.  Mtoivs  f.iivcrril  Media  (Academic, 

New  Yolk,  ItHiOl.  ('hap.  IV. 

5A.  Snmnieifeld,  Optics,  Vol.  IN*  of  t.eetures  ttn  Theoretical 
l*h\sics  (Academic.  Now  York. 

I).  Pergtand.  tFFE  Spectrum,  It  (I9f*{0, 

\i.  W.  Cooley  and  ,1.  W.  Tukey,  Math.  Cnmput,  19.  297  (UMiM. 

'  N.  ttienner.  IBM  Corp.  Contributed  Program  Library,  Code 
No.  .ISO,  M.’l.  l.ont  tunpublished). 

L.  K rm| niff,  Oeophysles  24.  .'10  (t9.r.9K 
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*1  tow  ever,  it  must  lie  pointer!  out  that,  while  this  is  true,  con¬ 
ventional  imaging  also  becomes  much  more  effective  as  the 
angular  separation  increases.  In  fact,  conventional  imaging 
improves  more  rapidly,  othe r  things  being  equal.  Thus,  the 
WSET  is  increasingly  preferable  in  comparison  with  other 
Methods  as  the  angular  separation  is  reduced. 
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EXPLANATION 


Younger  alluvium  and  colluvium 
Unconsolidated  gravel,  sand,  silt,  and  clay  in  main  stream 
beds  and  large  talus  deposits  on  slopes 


Qalo 


Older  alluvium  and  colluvium 
Chiefly  carbonate-cemented  or  unconsolidated  sand  and 
granule  and  pebble  conglomerate 

~Toev 

~~To» 

~Toeb 

Oak  Spring  formation 

Tosw,  gray  pumiceous  crystal  vitric  welded  tuff . 

Tos,  white,  red,  yellow,  pink-gray  pumiceous  nonwelded 
crystal  vitric  tuff,  zeolitized,  especially  in  lower 
part. 

Tosb,  widespread  basal  unit;  reddish-brown  tuffaceous 
rubble  breccia  and  conglomerate,  and  conglomeratic 
tuff 


!h 
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Granite,  light-gray,  equigranular,  fine-  to  coarse-grained; 
aplite  to  pegmatite;  occurs  as  dikes  (gd)  and  sills  (gs); 
"a"  prefixed  to  gd  or  gs  indicates  feldspathic,  argillic, 
sericitic,  or  silicic  alteration  prominent  in  dike  or 


TKg 


Granodiorite,  light-  to  greenish-gray,  predominantly  medium 
grained,  some  fine-grained;  occurs  as  a  stock  and  sills 


TKq 


Quartz  monzonite,  light-  to  medium-gray,  porphyritic,  fine- 
to  medium- grained;  occurs  as  a  stock 


Xf 


PPt 


Tippipah (?)  1 ime stone 

Dolomite,  medium-  to  dark-gray,  aphanitic  to  finely  crystal 
line,  thick-  to  very  thick  bedded,  nonfossiliferous.  Es¬ 
timated  minimum  thickness  200  feet 


Me 


Meu 

Mem 

Mel 


Eleanu  formation 

Me,  undifferentiated 

Meu,  argillite,  chiefly  very  dark  gray,  red-purple  or 
buff,  laminated  to  thin-bedded.  Limestone,  light-  to 
very  dark-gray,  finely  to  coarsely  crystalline,  partly 


PENNSYLVANIAN  CRETACEOUS  (? )  OR  TERTIARY  ( ? )  TERTIARY 


PENNSYLVANIAN  CRETACEOUS  (? )  OR  TERTIARY(?) 

AND  PERMIAN ( ? ) 
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PLANATION 
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Note:  Not  indicated  specifically  are  broad  areas  ol  the 
stock  that  have  undergone  pervasive  hydrothermal  py- 
rite  mineralization  and  argillic  alteration  of  plagio- 
clase  and(or)  chloritic  alteration  of  blotite  Also 
not  specifically  indicated  is  the  extensive  silica- 
tion  of  the  marble  along  fractures  near  the  edg' s 
of  the  stock  The  following  areas  of  alteration  have 
been  mapped  mainly  in  the  stock  except  as  noted 


tt 


sin 


Silicic  alteration,  strong(ss)  and  moderate  (sm).  Ex¬ 
cludes  individually  mappable  veins  or  groups  ol  veins 


'/ 


Quartz  vein  or  group  of  veins 


B  V 
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Sericitic 

Sericite  (white  mica)  is  predominant  or  prominent  mineral 


Hematitic 

lematite  in  prominent  amounts.  Generally  associated  with 
masses  of  altered  igneous  rock  in  which  quartz  is  prom¬ 
inent  or  along  siliceous  or  silicate  veins  in  the 
marble 
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/  V 


agd 


Feldspathic 

(without  letter  symbol  on  map) 

Lens-shaped  feldspathized  masses  of  igneous  rock,  commonly 
argillized  and  subordinately  sericitized,  silicified, 
and  pyritized.  The  masses  are  commonly  stained  red, 
purple,  yellow,  and(or)  brown  by  iron(?)  oxides 
Because  of  the  similarity  of  the  size  and  shape  at  map 
scale,  the  granite  dikes  affected  by  the  same  altera¬ 
tion  are  indicated  by  "a" 


Tactite 

Predominantly  garnet  with  subordinate  quartz,  epidote, 
green  mica(?),  limonite,  calcite,  sulfide,  tungstate, 
and  copper  carbonate  minerals 


Contact,  showing  dip.  Long  dashed  where  approximately 
located,  short  dashed  where  gradational  or  inferred, 
dotted  where  concealed 
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Geology  mapped  n  19^)9  and  1960 


uy±,  uoj.onu.-ce,  wmme  zo  medium- gray,  poorly  to 
moderately  foliated;  contains  some  nonpersistent 
interlayers  up  to  75  feet  thick  with  numerous 
nodules 


■Cs 

Stirling(?)  quartzite 

Quartzite,  light  to  medium  pinkish-gray  and  dark-gray, 
medium-grained,  poorly  bedded  with  rare  thin  beds 
of  dark-  to  medium- gray  phyllite 


In  order  to  facilitate  early  publication  of 
this  map,  some  of  the  normal  cartograpt  ic 
refinements  have  been  omitted 
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